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PREFACE 


It  is  now  well  recognized  that  certain  types  of  macromolecules , 
when  added  in  relatively  small  concentrations  to  solutions  flowing  in 
turbulent  motion,  possess  the  property  of  reducing  the  friction  drag 
of  the  fluid  in  which  they  are  suspended  or  dissolved.   As  a  result, 
the  addition  of  these  substances  can  significantly  increase  the  capacity 
of  a  system  to  transmit  fluids  at  normally  turbulent  flow  rates. 

The  feasibility  of  utilizing  this  effect  to  reduce  the  time 
required  to  perform  certain  critical  refueling  operations  encountered 
in  military  logistics  support  prompted  this  investigation.   That  the 
tactical  vulnerability  of  a  combatant  force  is  at  its  peak  during 
periods  of  logistic  replenishment  is  well  recognized  by  military  logis- 
ticians.   This  is  particularly  true  during  periods  involving  underway 
refueling  operations  from  a  mobile  replenishment  group  at  sea,  midair 
refueling  of  aircraft  and  the  refueling  of  aircraft  aboard  aircraft- 
carriers  at  sea. 

The  original  scope  of  this  investigation  envisioned  an  examina- 
tion of  the  drag  reducing  effects  of  a  number  of  viscoelastic  polymers 
in  several  of  the  more  common  military  fuels  in  use.   It  was  further 
recognized  that  the  effect  these  additives  would  have  on  the  quality  of 
the  fuel  would  be  of  paramount  interest.   However,  because  of  the  time 
required  to  design,  assemble  and  calibrate  a  flow  system  which  would  be 
suitable  for  investigation  of  the  drag  reduction  phenomena,  this  scope 
was  limited  to  the  examination  of  a  single  high  molecular  weight  poly- 
isobutylene  (Vistanex  L-200)  added  to  military  jet  fuel  (JP-4).   The 
effect  of  this  additive  on  fuel  quality  was  not  investigated. 

A  variable  pressure,  single  pass  flow  system  was  designed  and 
constructed  for  use  in  this  study.   This  type  of  system  was  selected  in 
order  to  maintain  similarity  with  tactical  refueling  operations. 
Further,  the  effects  of  mechanical  degradation  could  best  be  observed 
in  a  single  pass  system  since  the  fluid  makes  only  one  pass  through  the 
test  section.   Finally,  flow  behavior  over  a  wide  range  of  the  turbulent 
flow  regime  could  be  observed  yet  fluid  consumption  and  equipment  cost 
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could  be  minimized.   The  apparatus  was  found  to  perform  well  in  both 
the  laminar  and  turbulent  flow  regimes  when  tested  with  liquids  of  well 
known  fluid  properties. 

All  drag  reduction  experiments  were  accomplished  by  observing 
the  flow  behavior  of  solutions  containing  50  wppm  of  the  high  molecular 
weight  Vistanex  L-200  in  military  jet  fuel.   The  results  of  the  study 
indicate  that  a  marked  reduction  in  the  viscous  drag  of  jet  fuel  flow- 
ing in  the  turbulent  regime  can  be  achieved.   It  is  anticipated  that 
the  time  required  to  transmit  fuel  could  be  better  than  halved. 

The  polymer  tested  was  found  to  be  extremely  sensitive  to 
mechanical  degradation.   While  the  results  of  this  study  are  not  conclu- 
sive, a  striking  similarity  between  what  has  been  termed  "the  diameter 
effect"  and  the  effects  of  mechanical  degradation  was  observed.   A  cor- 
relation suitable  for  dynamic  scale-up  to  full  sized  flow  systems  was 
frustrated  by  this  degradation/diameter  effect.   Further,  it  was  found 
that  mechanical  devices  producing  extremely  high  localized  shear  rates 
could  completely  destroy  the  drag  reducing  effectiveness  of  the  additive , 
Therefore,  it  was  concluded  that  this  type  of  a  device  should  be  avoided 
in  any  full  scale  flow  system  or,  alternatively,  injection  of  the 
additive  should  be  accomplished  downstream  from  such  a  device. 

There  were  strong  indications  that  a  suppression  of  turbulent 
fluctuations  might  well  be  the  main  characteristic  of  drag  reduction. 
Additional  work  with  more  precise  equipment  was  recommended  which  would 
serve  to  verify  this  tentative  conclusion. 
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CHAPTER  I 
INTRODUCTION 

The  objective  of  this  study  was  to  investigate  the  effects  on 
fluid  flow  caused  by  the  addition  of  a  minute  quantity  of  a  high 
molecular  weight  polymer  to  military  jet  fuel  (JPH).   The  flow  was 

observed  in  the  turbulent  flow  regime  corresponding  to  a  range  of 

4      5  .        . 

Reynolds  numbers  from  10   to  10  .   In  particular  this  thesis  is  con- 
cerned with  the  characteristics  of  frictional  drag  (or  its  reduction)  as 
influenced  by  such  an  additive  under  conditions  of  flow  in  cylindrical 
tubes.   The  primary  emphasis  of  this  study  has  been  to  seek  the  factors 
which  tend  to  influence  the  effectiveness  of  such  additives  in  a  common 
fuel  in  actual  use,  rather  than  to  search  for  the  mechanisms  which  might 
be  responsible  for  the  cause  of  drag  reduction.   Further,  it  was  recog'-. 
nized  that  time  would  permit  only  a  very  limited  examination  of  the 
influence  of  these  additives  since  it  would  first  be  necessary  to  design, 
assemble  and  calibrate  a  flow  system  suitable  for  the  investigation  of 
the  drag  reduction  phenomenon.   However,  it  was  hoped  that  this  effort 
might  become  the  basis  for  future  work  in  this  area  since  a  calibrated 
system  would  be  readily  available. 

LITERATURE  SURVEY  AND  DISCUSSION 

(a)   Drag  Reduction  Defined   The  term  "drag  reduction"  or 
"fluid  friction  reduction"  has  been  used  to  describe  a  rather  broad 
spectrum  of  rheological  phenomena  since  first  introduced  in  the  technical 
literature  in  1949.   It  is  therefore  necessary  to  consider  two  quite 
distinct  classifications  of  fluids,  purely  viscous  and  viscoelastic,  with 
which  the  term  is  frequently  associated  in  order  to  clarify  the  defini- 
tion as  used  in  this  thesis.   The  first  classification  consisting  of 
fluids  commonly  described  as  "purely  viscous"  can  be  further  subdivided 
into  two  categories.   The  fluids  of  the  first  category  are  said  to  obey 
Newton's  Law  of  Viscosity  and  are  broadly  classified  as  "Newtonian 
Fluids."   The  behavior  of ' such  fluids  can  be  described  by  the  general 
power  law  or  "Ostwald-de  Waele"  model  (3): 
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T  =  -  m  (dv  /dy)n  [1-1] 

where  n  =  1  and  m  is  a  constant  and  m  =  y 

This  relationship  states  that  the  shear  force  per  unit  area  (shear 
stress)  varies  directly  with  the  negative  of  the  velocity  gradient.   A 
more  useful  relationship,  since  it  describes  the  shear  stress/shear 
rate  characteristics  of  laminar  flow  in  cylindrical  tubes  in  terms  of 
measureable  or  derived  quantities,  is: 

x      =   DAP/4L  =  y  (8V/D)  [1-2] 

w 

or 

RAP/2L  =  y  (4Q/ttR3)  [1-3] 

Results  of  extensive  experimental  tests  of  these  fluids  in  turbulent 
flow  through  sand-roughed  pipes  were  presented  by  Nikuradse  in  1933. 
Friction  factor  data  are  conventionally  presented  as  functions  of 
Reynolds  numbers  in  dimensionless  plots  using  logarithmic  scales  with 
relative  roughness  as  a  secondary  variable.   Perhaps  the  most  well 
known  plots  of  this  type  are  the  "Moody  Diagrams"  (14)  found  in  common 
engineering  handbooks.  (9) 

Fluids  in  the  second  category,  "purely  viscous,"  are  classified 
as  "non-Newtonian"  and  can  generally  be  described  by  the  power  law  rela- 
tion, equation  [1-1],  but  without  the  restrictions  that  n  =  1  and  m  be 
constant.   Similarly,  in  terms  of  physical  quantities  more  amenable  to 
direct  measurement: 

T   =  DAP/4L  =  m  [2(3  +  1/n)  V/D]n  [1-4] 

w 

After  measuring  AP  as  a  function  of  V  (or  flow  rate  from  which  an  aver- 
age velocity  can  be  derived)  a  plot  of  the  variables  DAP/4L  versus 
8V/D  on  log-log  paper  can  be  constructed  which  will  reveal  to  which  of 
these  two  viscous  categories  a  fluid  belongs.   The  popularity  of  the 
Ostwald-de  Waele  power  law  as  a  means  of  describing  viscous  fluid 
behavior  now  becomes  readily  apparent  since  it  has  been  shown  to  cor- 
rectly portray  the  relationship  between  shear  stress  and  shear  rate  for 
many  fluids  over  the  broad  range  of  shear  rates  commonly  encountered. 
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Extensions  of  the  experiments  conducted  by  Nikuradse  which 
include  the  behavior  of  non-Newtonian  fluids  have  been  the  subject  of 
more  recent  investigators.   (4,  22)   By  using'  a  "generalized"  Reynolds 
number  which  has  been  shown  to  be  a  function  of  power  law  coefficients, 
m  and  n,  Dodge  and  Metzner  (5)  have  revealed  the  nature  of  the  friction 
factor,  "generalized"  Reynolds  number  relationship  for  flow  of  non- 
Newtonian  fluids  through  smooth  cylindrical  pipes.   The  results  of  the 
Dodge  and  Metzner  studies  indicate  that : 

(a)  The  generalized  Reynolds  number  at  which  transition 
from  laminar  flow  to  turbulent  flow  is  observed  increases  slowly  from  a 
value  of  2100  for  Newtonian  fluids  as  the  flow  behavior  index,  n, 
decreases.   (n  is  customarily  defined  as  the  index  of  non-Newtonian 
behavior. ) 

(b)  At  any  given  generalized  Reynolds  number  in  the  turbu- 
lent flow  regime  there  will  be  a  decrease  in  the  friction  factor 
observed  corresponding  to  a  decrease  in  the  flow  behavior  index,  n. 

These  results  might  suggest  behavior  not  totally  unlike  that  associated 
with  the  drag  reduction  phenomenon. 

Since  the  friction  factors  for  purely  viscous  non-Newtonian 
fluids  usually  will  be  lower  than  those  of  Newtonian  materials  if  n<l,  it 
would  appear  as  if  the  additions  of  solids  or  polymeric  materials  to  a 
Newtonian  solvent  to  yield  a  non-Newtonian  fluid  might  serve  to  reduce 
pressure  drop  under  turbulent  conditions.   In  fact,  however,  the  thicken- 
ing action  promoted  by  these  additives  generally  will  outweigh  the  fric- 
tion factor  reduction  so  that  the  actual  pressure  drop  required  to 
develop  a  certain  flow  rate  will  almost  always  be  increased.   As  correctly 
commented  upon  by  Savins  (21)  and  Fabula  and  Hoyt  (6),  it  is  important 
to  recognize  the  difference  between  non-Newtonian  behavior  and  the  drag 
reduction  phenomenon,  which  is  to  be  discussed  shortly. 

The  high  molecular  weight,  random-coiling  polymers  in  solution 
appear  to  exhibit  elastic  characteristics  in  addition. to  the  viscous 
characteristics  which  describe  the  behavior  of  most  fluids.   Hence,  the 
second  broad  classification  of  fluids  to  be  considered  is  termed 
"viscoelastic."   It  might  be  well  to  recall  at  this  point  that  the  prop- 
erties of  elasticity  and  viscosity  are  defined  in  terms  of  the  response 
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of  a  fluid  to  an  applied  stress.   The  elastic  response  describes  the 
tendency  of  a  fluid  to  recover  its  former  configuration  upon  removal 
of  an  applied  stress.   This  recovery  would  be  complete  were  it  not  for 
the  viscous  flow  which  begins  when  the  stress  is  applied.   The  viscous 
flow  is  never  reversed  and  may  continue  once  the  stress  has  been  removed 
owing  to  the  inertial  properties  associated  with  the  fluid.   Thus,  an 
essential  distinction  between  elasticity  and  viscosity  is  that  the 
former  describes  a  reversible  process  while  the  latter  describes  an 
irreversible  property. 

As  might  well  be  expected,  very  dilute  solutions  containing  these 
viscoelastic  polymers  exhibit  characteristics  very  closely  identified 
with  the  fluid  properties  described  as  "Newtonian"  by  conventional  vis- 
cometry  i.e.  ,  n  =  1,  m  =  y  in  Eq.  [1-1].   Owing  to  the  very  low  concen- 
trations of  polymer  additives,  the  measured  density  and  viscosity  of  a 
solution  are  very  nearly  the  same  as  those  of  the  pure  solvent.   At- 
higher  levels  of  concentration,  they  appear  to  impart  non-Newtonian 
characteristics  and  while  the  assignment  of  "viscosity"  as  a  fluid 
property  to  such  fluids  is  somewhat  nebulous ,  a  thickening  of  the  fluid 
is  readily  apparent.   While  investigating  the  characteristics  of  non- 
Newtonian  fluids,  both  Dodge  and  Metzner  (M-,  5)  and  Shaver  and 
Merrill  (22)  selected  some  additives  which  are  now  known  to  produce 
viscoelastic  solutions  .   An  apparent  anomalous  behavior  was  observed  for 
the  viscoelastic  fluids  when  compared  to  the  purely  viscous  fluids  with 
approximately  the  same  index  of  non-Newtonian  behavior. 

It  is  at  this  point  that  authorities  appear  to  disagree  on  the 
meaning  of  "drag  reduction."   Shin  (23)  contends  that  if  a  polymer  solu- 
tion is  "concentrated,"  it  will  display  considerable  gradient  dependency 
of  viscosity  (the  characteristic  most  commonly  associated  with  non- 
Newtonian  behavior)  and  that  the  assignment  of  viscosity  for  purposes 
of  evaluating  a  Reynolds  number  is  always  equivocal.   This  suggests  that 
the  "flow  anomaly"  associated  with  the  concentrated  viscoelastic  solu- 
tions observed  by  Shaver,  Dodge,  et  al.,  might  well  stem  from  an  imper- 
fect model  used  to  evaluate  viscosity  and  generalized  Reynolds  number 
rather  than  an  actual  reduction  in  friction  stress.   In  contrast,  the 
Newtonian  characteristics  of  the  solvent  are  retained  when  dilute  solu- 
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tions  containing  viscoelastic  polymers  are  formed  and  the  classical 
properties  associated  with  viscosity  and  Reynolds  number  can  be  unequiv- 
ocally evaluated  in  a  conventional  sense.   Further,  in  very  dilute 
polymeric  solutions,  the  fluid  properties  generally  used  to  describe 
fluid  flow  are  very  nearly  equal  to  the  properties  of  the  pure  solvent 
when  subjected  to  conventional  measurement  techniques.   Thus,  the  depar- 
ture from  conventional  flow  behavior  associated  with  these  extremely 
dilute  solutions  can  very  clearly  be  labeled  "anomalous."   This  is  not 
to  suggest  that  these  same  characteristics  do  not  occur  in  concentrated 
viscoelastic  solutions.   Indeed,  the  observations  of  Dodge  (4)  and 
Shaver  and  Merrill  (22)  and  more  recently  those  of  Savins  (21)  and 
Rodriguez,  Zakin  and  Patterson  (20)  strongly  suggest  that  frictional 
stress  reducing  characteristics  separate  and  apart  from  non-Newtonian 
behavior  of  purely  viscous  fluids  are  present  in  these  concentrated 
solutions.   However,  from  a  practical  standpoint,  the  final  criterion 
of  the  effectiveness  of  such  an  additive  must  be  measured  by  its  fric- 
tion reducing  abilities  as  compared  to  the  drag  of  the  pure  solvent  in 
which  it  is  dissolved.   Thus,  whether  a  final  solution  is  considered 
Newtonian  or  non-Newtonian,  dilute  or  concentrated,  it  is  the  frictional 
drag  of  that  solution  as  compared  to  the  frictional  drag  of  the  pure 
solvent  which  shall  be  used  as  the  criterion  to  define  drag  reduction 
in  this  thesis. 

(b)  Causes  and  Effects  of  Drag  Reduction  Phenomena   While 
there  have  been  a  number  of  attempts  to  isolate  the  precise  mechanisms 
which  give  rise  to  the  drag  reduction  phenomena  noted  in  the  literature, 
the  exact  mechanism  is  still  not  known.   This,  of  course,  is  understand- 
able considering  that  a  microscopic  phenomenon  is  being  interpreted  on 
the  basis  of  macroscopic  data  fitted  to  a  macroscopic  or  microscopic 
model. 

The  questions  to  be  examined  if  one  is  to  gain  insight  into  the 
mechanisms  which  give  rise  to  drag  reduction,  can  now  be  stated  more 
explicitly  as :   how  the  macroscopic  fluid  properties  and  flow  behavior 
are  related  to  the  microscopic  configuration  and  motion  of  the  molecules, 
and,  how  the  measurable  macroscopic  quantities  such  as  fluid  density, 
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viscosity  and  flow  velocity  can  be  shown  to  be  consistent  with  a  macro- 
scopic or  microscopic  model  of  fluid  flow. 

The  earliest  published  technical  report  available  which  attributes 
drag  reduction  to  dilute  polymer  solutions  appears  to  be  that  of  B.A. 
Toms  appearing  in  the  Proceedings  of  the  First  International  Rheological 
Congress.   (25)   Indeed,  because  of  the  pioneering  research  he  conducted 
in  this  field,  it  is  often  referred  to  as  the  "Toms  Effect."   (7)   In  his 
paper,  Toms  presented  data  which  indicated  an  observed  increase  in  the 
rate  of  flow  of  monochlorbenzene  at  a  constant  pressure  gradient  caused 
by  the  addition  of  polymethylmethacrylate.   The  remarkable  feature  of 
these  results  was  that  this  increased  flow  rate  was  a  phenomenon 
observed  only  in  the  nonlaminar  flow  region. 

Oldroyd  (15)  concurrently  presented  a  paper  suggesting  that  these 
unexpectedly  high  flow  rates  observed  by  Toms  be  explained  in  terms  of 
anomalous  behavior  of  a  thin  laminar  layer  at  the  tube  wall  when  the 
mainstream  of  flow  becomes  turbulent.   Oldroyd  reasoned  that,  if  the 
fluid  is  a  solution  containing  a  long  chain  linear  polymer,  an  external 
constraint  may  be  imposed  by  the  wall  of  the  flow  system  on  the  ways  in 
which  these  long  chain  molecules  can  rotate.   It  might  therefore  be 
possible  that  an  abnormally  mobile  pseudo-laminar  sublayer,  or  boundary, 
of  a  thickness  comparable  with  molecular  dimensions  could  exist  at  the 
wall.   This  effect  would  then  be  macroscopically  evident  as  an  effective 
slippage  at  the  wall.   The  flow  rates  as  determined  for  conventional 
Newtonian  fluids  would  have  to  be  increased  to  account  for  this  slippage. 

Although  the  proposed  velocity  of  slip  might  well  explain  the 
macroscopic  effect   of  drag  reduction,  a  departure  from  the  accepted 
model  of  fluid  flow  is  implied  if  this  theory  is  used  to  describe  the 
microscopic  cause.      In  contrast  to  the  model  used  to  describe  ordinary 
flow  of  a  Newtonian  fluid,  it  was  visualized  that  the  long  chain  poly- 
mers moved  with  translational  motion,  hence  would  slide  or  slip  along 
the  side  of  the  wall.   Subsequent  investigation  however,  has  failed  to 
reveal  any  translational  slippage.   To  the  contrary,  flow  experiments 
conducted  in  both  the  laminar  and  turbulent  flow  regimes  appear  to 
confirm  a  flow  model  based  on  a  zero  velocity  at  the  wall.   Despite  the 
fact  that  Oldroyd' s  slip  mechanism  as  a  cause  for  drag  reduction  has  not 


been  given  serious  consideration  for  some  time,  there  are  relatively 
current  texts  which  continue  to  refer  to  slip  as  being  "important  in 
some  non-Newtonian  flow  problems."   (3) 

A  second  concept  which  might  explain  the  mechanism  or  cause  of 
drag  reduction  was  first  advanced  by  Davies,  Ward,  Atkinson  and  Blair 
(17)  in  reaction  to  Oldroyd's  boundary  layer  theory.   These  gentlemen 
suggested  the  possibility  that  the  Toms  Effect  is  tied  to  the  elasticity 
of  polymer  solutions,  and  although  variations  of  this  theory  appear  to 
be  the  most  widely  accepted  today,  the  exact  mechanism  by  which  the 
elastic  properties  influence  drag  reduction  is  still  not  fully  under- 
stood.  In  general,  this  theory  holds  that  the  elastic  properties  of 
the  fluid  absorb  the  kinetic  energy  of  turbulent  motion,  thereby  reduc- 
ing the  fluid  transmission  energy  lost  to  turbulence.   In  other  words, 
the  irreversible  friction  losses  due  to  turbulence  are  reduced  because 
the  energy  associated  with  turbulent  fluctuations  has  been  reversibly 
stored  in  the  fluid. 

A  third  mechanism,  the  formation  of  a  viscosity  gradient,  was 
apparently  originally  proposed  by  Shaver  and  Merrill  (22)  though  it  is 
discussed  in  some  detail  by  Shin  (25)  under  the  heading  "Anisotropic 
Viscosity."   This  theory  stems  from  Merrill's  suggestion  that  the  local 
viscosity  normal  to  the  direction  of  shear  might  well  be  substantially 
increased  through  elongation  of  the  randomly  coiled  macromolecules  even 
though  the  viscosity  parallel  to  the  direction  of  shear  is  very  nearly 
identical  to  that  of  the  solvent.   Thus,  the  high  perpendicular  vis- 
cosity can  reduce  turbulent  fluctuations  in  that  direction.   Since  the 
turbulent  fluctuations  are  randomly  directed,  the  orientation  of  the 
elongated  molecules  will  also  be  random. 

During  the  early  1960 's,  the  Navy  became  particularly  interested 
in  the  drag  reduction  phenomenon  because  of  its  possible  influence  on 
ship  and  torpedo  performance.   Fabula  and  Hoyt  (6)  of  the  Navy's 
Ordnance  Test  Station  at  China  Lake,  California,  were  particularly 
instrumental  in  surveying  many  water  soluble  drag  reducing  additives. 
As  a  result  of  their  experiments  conducted  with  aqueous  polymer  solu- 
tions, general  rules  were  developed  as  to  the  type  of  material  likely 
to  be  effective.   In  general,  it  was  concluded  that  at  least  three 
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significant  parameters  affect  the  ability  of  a  polymer  to  lower  the 
turbulent  frictional  resistance  of  the  fluid  in  which  it  is  dissolved: 
linearity,  molecular  weight,  and  solubility.   The  "long  chain"  mate- 
rials have  an  essentially  unbranched  molecular  structure  of  high 
molecular  weight  when  added  to  solvents  in  which  they  can  be  fully 
dissolved,  appear  to  produce  the  most  marked  drag  reduction  effect. 
A  review  of  the  work  of  Fabula  and  Hoyt  also  reveals  that:   (a)  drag 
reduction  effectiveness  tends  to  improve  with  increased  concentrations 
up  to  a  certain  concentration  level,  after  which  any  further  increase 
in  concentration  level  will  result  in  decreased  effectiveness;   (b)  at 
concentration  levels  less  than  the  "optimum"  concentration  for  any 
given  additive,  the  friction  factor  versus  Reynolds  number 
curves  are  roughly  parallel  to,  but  lower   than,  the  friction  factor 
curve  for  the  pure  solvent  flowing  in  the  turbulent  flow  regime ; 
(c)  at  concentration  levels  above  the  optimum  concentration  the  friction 
factor  curves  are  roughly  parallel  to,  but  higher-   than,  the  extended 
laminar  fluid  flow  curve  of  the  pure  Newtonian  solvent.   Additional 
effects  which  appear  to  be  more  related  to  the  flow  system  than  to  the 
fluid  are  also  commonly  observed: 

(b)   For  a  given  solution  which  exhibits  drag  reduction 
characteristics  the  drag  reduction  effect  is  more  pronounced  as  the 
diameter  of  the  flow  system  is  increased. 

(b)   For  a  given  solution  and  flow  system,  the  drag  reduc- 
tion effect  becomes  more  pronounced  as  flow  rates  are  increased. 
However,  there  appears  to  be  a  maximum  tolerable  shear  rate  associated 
with  each  flow  system  beyond  which  the  drag  reducing  effectiveness 
appears  to  decline. 

The  more  recent  work  of  Shin  (25)  has  served  to  quantify  some  of 
these  qualitative  observations,  particularly  with  respect  to  defining 
the  significant  characteristics  of  a  "dilute,"  "concentrated,"  "critical" 
and  "optimum"  solution.   Shin's  results  suggest  that  the  key  to  the 
question  of  whether  the  polymer  solution  is  dilute  or  concentrated 
depends  on  the  critical  concentration  which  corresponds  to  the  condition 
that  each  polymer  molecule  (having  been  assigned  an  effective  spherical 
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diameter  which  is  of  the  order  of  magnitude  of  its  r.m.s.  end-to-end 
distance)  is  brought  to  a  spherical  packing  of  74%.   In  other  words, 
the  random  coils  touch  each  other,  but  are  otherwise  unchanged.   When 
the  polymer  concentration  is  less  than  the  critical  concentration,  it 
is  considered  "dilute."   It  is  when  a  solution  is  "concentrated"  that 
it  will  display  considerable  gradient  dependency  of  viscosity  or  non- 
Newtonian  behavior.   On  the  other  hand,  polymer  solutions  less  than  one 
third  of  the  critical  concentration  appear  to  display  either  negligible 
or  undetectable  gradient  dependency  of  viscosity.   While  there  was  no 
direct  correlation  achieved  between  the  optimum  concentration  and 
critical  concentration,  the  optimum  concentration  was  observed  to  be 
generally  one  tenth  to  one  hundredth  of  the  critical  concentration. 
For  example,  the  critical  concentration  of  polyisobutylene  L-200  in 
cyclohexane  was  found  to  be  520  wppm;  the  optimum  drag  reduction  effect- 
iveness over  the  limited  shear  rates  investigated  was  found  to  be 
35  wppm.   As  was  emphasized  by  Shin,  it  is  not  commonly  recognized  that 
the  critical  concentration  of  solutions  of  polymers  having  molecular 
weights  in  the  millions  is  small.   Indeed,  0.1%  polymer  solutions  are 
usually  "concentrated"  as  defined  by  Shin. 

In  even  more  recent  work,  Rodriguez  et  al.  (20)  appear  to  have 
achieved  a  very  useful  correlation  between  drag  reducing  characteris- 
tics for  turbulent  flow  in  a  pipe  and  measurable  properties  of  the 
polymer  solution  and  flow  system.   They  have  successfully  correlated 
friction  reduction  behavior  of  a  viscoelastic  solution  with  a  modified 
"Deborah  number"  defined  as: 

N    -  V  t,/D°-2  [1-5] 

Db      1 

Fabula  et  al.  (7)  and  more  recently  Hershey  (8),  Astarita  (2) 
and  others  have  suggested  that  fluid  friction  reduction  occurs  when  the 
relaxation  time  of  the  polymer  molecules  in  solution  exceeds  a  character- 
istic flow  time  of  the  experiment.   The  relaxation  time  of  the  polymer 
solution  represents  the  relative  amounts  of  viscous  and  elastic  response. 
Characteristic  flow  time  has  been  taken  as  the  reciprocal  of  the  shear 
rate  at  the  wall.   The  problem  has  been  that  of  identifying  specific 
physical  quantities  susceptible  to  direct  measurement,  which  would  be 
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characteristic  of  fluid  relaxation  time.   Park  (13)  has  shown  that 
intrinsic  viscosity  will  not  yield  a  favorable  correlation  when  evalu- 
ating concentrated  polymeric  solutions.   Yet  the  basis  for  determining 
a  measure  of  relaxation  time  (i-,  )  in  the  Rodriguez  correlation  would 
appear  to  be  intrinsic  viscosity  (though  fluid  velocity  and  flew  system 
diameter  are,  of  course,  reflected  in  the  modified  Deborah  number). 
There  would  appear  to  be  some  disparity  between  the  results  of  Park  and 
Rodriguez.   It  would  be  interesting  to  investigate  the  drag  reducing 
capabilities  of  Carbopol  using  the  Rodriguez  correlation  since  this 
additive  would  appear  to  be  ideally  suited  to  promote  drag  reduction, 
though  none  has  been  observed  in  Carbopol  solutions  to  date! 

SCOPE  AND  LIMITATIONS  OF  EXPERIMENTAL  WORK 

The  original  scope  of  this  investigation  envisioned  an  examina- 
tion of  the  drag  reducing  effects  of  a  number  of  polymers  in  several  of 
the  more  commonly  used  military  fuels.   It  was  further  recognized  that 
the  effect  of  these  additives  would  have  on  the  quality  of  the  fuel 
would  be  of  paramount  interest.   However,  it  would  first  be  necessary  to 

design,  assemble  and  calibrate  a  flow  system  which  would  be  suitable 
for  investigation  of  the  drag  reduction  phenomena.   Because  of  the  time 
involved  in  this  phase  of  the  work,  the  original  scope  was  limited  to 
a  brief  examination  of  a  high  molecular  weight  additive  in  military  jet 
fuel  at  a  concentration  of  50  weight  parts  per  million. 

Time  also  placed  important  limitations  on  the  type  of  flow  system 
which  could  be  assembled  and  used  for  drag  reduction  experiments. 
Although  a  closed  loop,  recirculating  type  of  flow  system,  patterned 
after  those  of  Dodge  (4)  or  Melton  and  Malone  (12)  was  initially 
considered,  it  became  clear  that  the  time  required  to  obtain  and  assemble 
the  equipment  needed  for  such  an  elaborate  system  would  be  excessive. 
It  might  also  be  pertinent  to  note  that  systems  which  provide  for  a 
recirculation  of  the  polymeric  solution  will  not  permit  as  detailed  an 
examination  of  such  important  factors  as  mechanical  degradation  as  is 
permitted  in  a  single  pass  flow  system.   This  is  particularly  true  if 
recirculation  is  accomplished  by  a  mechanical  device  generating  high 
localized  shear  rates.   For  these  reasons,  it  appeared  that  a  single 
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pass  system  patterned  after  a  variable  pressure  capillary  tube  visco- 
meter (such  as  the  one  used  by  Dodge  for  viscometric  measurements) 
would  be  the  most  practical  system  which  could  be  assembled  in  a  reason- 
able time  period.   Hypodermic  tubes  for  use  as  test  sections  were  readily 
available  as  was  a  supply  of  bottled  nitrogen  gas  for  use  as  the 
driving  mechanism.   Since  it  was  desirable  to  design  a  flow  system 
capable  of  providing  a  means  of  observing  the  flow  behavior  in  the  well 
developed  turbulent  regime,  a  computerized  model  of  a  capillary  tube 
viscometer  was  developed  through  which  the  design  parameters  affecting 
fluid  flow  over  a  range  of  Reynolds  numbers  up  to  10   could  be  examined. 
The  maximum  driving  pressure,  of  course,  could  not  exceed  the  pressure 
provided  by  the  bottled  nitrogen.   Further,  the  largest  fluid  reservoir 
available  capable  of  withstanding  the  pressure  from  the  bottled  nitro- 
gen appeared  to  be  an  empty  nitrogen  cylinder  which  had  an  estimated 
capacity  of  between  10  and  15  gallons. 


CHAPTER  II 
DESIGN,  CONSTRUCTION,  CALIBRATION,  AND  OPERATION  OE  FLOW  SYSTEM 

PRELIMINARY  DESIGN  CONSIDERATIONS 

Previous  experimental  work  involving  observations  on  drag  reduc- 
tion have  generally  been  performed  on  three  distinctly  different  types 
of  flow  systems : 

(a)  The  immersed  rotating  disk  or  cylinder. 

(b)  The  closed  loop,  recirculating  flow  system. 

(c)  The  open  loop,  single  pass  flow  system  through  small 
capillary  tubes. 

Generally,  these  involved  a  type  of  rotational  shearing  boundary 
or  cylindrical  tube  flow  system,  each  of  which  has  its  own  distinct 
advantages  and  disadvantages.  This   study  was  undertaken  with  a  specific 
application  in  mind--the  feasibility  of  reducing  the  time  required  to 
perform  certain  critical  refueling  operations  encountered  in  the  logis- 
tics support  of  tactical  military  maneuvers.   A  characteristic  of  these 
operations  is  the  "one  shot"  transmission  of  fluid  through  relatively 
short  fuel  lines.   A  modification  of  the  latter  system  was  selected  for 
the  experimental  work  in  this  study  in  order  to  achieve  a  degree  of 
similarity  to  the  characteristics  of  tactical  military  logistical 
systems. 

The  system  consisted  of  a  set  of  stainless  steel  hypodermic  tubes 
which  could  be  connected  to  a  high  pressure  reservoir  into  which  the 
fluid  to  be  examined  was  loaded.   The  fluid  was  forced  out  of  the  reser- 
voir by  a  siphoning  technique  which  was  achieved  by  connecting  the 
reservoir  to  a  source  of  bottled  nitrogen  (see  Figure  2-1). 

An  energy  balance  on  a  capillary  tube  flow  system  reduces  to  the 
following  form,  considering  steady  state,  isothermal  flow  and  neglecting 
the  effects  of  hydrostatic  pressure: 

(P2-P1)  /p  +  (K/2gc)(V22  -  Vl2)  +  2fV22L/gcD  =  0      [2-1] 

where  points  1  and  2  refer  to  the  entrance  and  exit  to  the  flow  tube. 
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In  capillary  tube  measurements,  it  is  customary  to  measure  the  total 
pressure  applied  to  the  reservoir  chamber  (P  )  and  to  consider  the 
velocity  in  the  reservoir  (V  )  essentially  as  zero.   Equation  [2-1] 
then  reduces  to: 

AP  =  (pV2/2g  )  (K  +  4fL/D)  [2-2] 

The  applied  pressure  then  Is  used  to  overcome  viscous  resistance  to  flow 
within  the  tube  and  to  provide  the  kinetic  energy  of  the  stream  as  well 
as  the  entrance  and  exit  losses.   These  latter  forms  of  energy  are 
usually  accounted  for  In  the  variable  K.   For  ordinary  laminar  Newtonian 
flow  probably  the  most  widely  accepted  value  of  K  is  2.24  (4).   Unfor- 
tunately, the  appropriate  value  of  the  total  correction  to  account  for 
kinetic  energy  and  entrance  and  exit  effects  Is  somewhat  In  doubt  for 
non-Newtonian  and  viscoelastic  fluids  (24).   One  can  of  course,  place 
a  lower  limit  on  the  value,  that  of  1.0  for  a  kinetic  energy  correction 

corresponding  to  a  perfectly  flat  velocity  profile  (thus,  a  correction  of 

2 

pV  /2g  )  Dodge  (4)  arbitrarily  assigned  a  correction  corresponding  for 

°  .  ...  2 

the  combined  entrance,  exit  and  kinetic  energy  effects  of  1.5  pV  /g 

corresponding  to  a  value  of  K  equal  to  3. 

The  Ideal  solution  to  this  dilemma  Is  to  make  the  unknown  value 
of  K  insignificant  when  compared  to  the  4fL/D  term.   In  the  laminar 
regime — for  example  at  a  Reynolds  number  of  800 — f  takes  on  a  vlaue  of 
0.02.   A  length  to  diameter  ratio  of  1:1000  will  result  in  a  value  of 
4fL/D  =80.   Under  these  conditions,  the  widest  conceivable  error  for  K 
is  probably  less  that  1%.   If  kinetic  energy,  exit  and  entrance  effects 
are  completely  disregarded,  the  error  would  be  less  than  3%.   On  the 
other  hand,  In  the  turbulent  regime  with  a  drag  reduction  mechanism 
operative,  the  observed  Fanning  friction  factors  will  frequently  be  less 
than  0.002.   The  same  uncertainty  in  K  as  in  the  previous  example  will 
now  result  in  an  error  of  almost  8%.   Further,  It  was  found  that  an  L/D 
ratio  in  excess  of  200  could  not  be  tolerated  if  the  goal  of  observa- 
tions extending  to  a  Reynolds  number  of  10  was  to  be  achieved  for  the 
full  class  of  fuels  originally  envisioned.   The  uncertainty  in  K  intro- 
duced by  this  new  criterion  could  perpetuate  an  error  of  as  much  as  20%. 
Specially  designed  flow  tubes  with  pressure  ports  located  well  downstream 
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from  the  entrance  to  the  tube  obviated  the  need  to  consider  these 
corrections . 

It  might  be  pertinent  to  note  that  drag  reduction  experiments 
performed  in  a  capillary  tube  rheometer  at  the  University  of  Minnesota 
were  reported  in  1964  (19).   The  flow  tubes  were  1000  diameters  in 
length  in  all  cases.   The  report  indicates  that  this  "arbitrary  selec- 
tion of  length  was  based  on  the  recognition  that  fluid  characterization 
errors  contributed  by  entrance,  exit,  etc.,  diminished  with  increasing 
length  while  errors  contributed  by  shear  degradation  might  increase  with 
length."  However,  the  highest  Reynolds  number  achieved  in  this  apparatus 
with  drag  reduction  effects  operative  was  approximately  25,000. 

A  computerized  model  of  a  capillary  tube  flow  system  was  developed 
through  which  preliminary  design  concepts  could  be  tested  for  feasibility 
and  practicality. 

The  computer  model  used  for  preliminary  design  was  expanded  as 
the  final  design  of  the  actual  flow  system  progressed,  so  that  the  effect 
of  each  design  feature  could  be  tested  in  advance  of  final  construction. 
Therefore,  computer  program  number  1  in  Appendix  A  accounts  for  the 
design  characteristics  peculiar  to  the  permanently  installed  portions  of 
the  flow  system.   This  program  can  be  used  for  two  different  purposes: 

(a)  To  provide  the  reservoir  pressure  which  will  be  required 
to  operate  the  flow  system  at  a  given  Reynolds  number. 

(b)  To  test  the  effect  that  variations  of  flow  tube  char- 
acteristics such  as  flow  tube  diameter,  entrance  length,  exit  length,  and 
test  section  length  will  have  on  the  performance  of  the  flow  system. 

The  former  feature  is  designed  to  assist  the  operator  of  the  flow 
system  by  providing  the  minimum  weight  or  volume  of  fluid  required  to 
meet  the  criterion  that  the  flow  time  for  the  run  last  at  least  100 
seconds,  and  by  providing  the  pressure  regulator  setting  required  to 
achieve  the  flow  regime  (at  a  given  Reynolds  number)  desired.   The  latter 
feature  is  provided  to  serve  as  a  guide  to  the  design  and  expected  per- 
formance of  new  flow  tubes  to  be  used  with  this  flow  system.   In  either 
case,  if  the  system  limitations  of  available  pressure  or  reservoir  capac- 
ity are  exceeded,  the  program  will  automatically  signal  such  a  condition 
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to  the  operator.   A  sample  write  out  from  this  program  is  included  in 
Table  2-1. 

DISCUSSION  OF  SUPPORT  EQUIPMENT  USED  AND  CALIBRATION  TECHNIQUES  EMPLOYED 

Fluid  flow  rates  were  determined  by  timing  the  flow  of  a  pre- 
determined quantity  of  fluid  with  a  Galco  stop  watch  (Chemical  Engineer- 
ing Department  #110).   The  watch  dial  indicated  time  in  increments  of 
0.2  seconds  but  there  was  a  sufficient  spread  between  marked  increments 
to  allow  interpolation  to  the  nearest  0.1  second.   A  predetermined 
criterion  was  imposed:  that  a  run  must  take  at  least  100  seconds  to 
insure  that  the  precision  of  the  time  measurement  was  1  part  in  1000. 
The  decision  as  to  whether  a  volumetric  or  weight  measurement  would  be 
taken  using  100  and  500  ml  graduated  cylinders,  a  2000  ml  volumetric 
flask  and  a  4000  Erlenmeyer  flask  which  was  calibrated  with  the  2000  ml 
flask  and  marked  at  the  4000  ml  level.   When  it  was  determined  that  the 
rate  of  flow  would  be  too  great  to  permit  a  time  measurement  of  4000  ml 
in  at  least  100  seconds,  a  weight  measurement  of  fluid  flow  was  taken 
using  a  Fairbanks  Morse  and  Company  200  pound  capacity  scale,  serial 
number  G  403886,  code  528.   This  scale  contained  marked  increments  of 
1/16  pound  with  a  sufficient  spread  between  the  marked  increments  to 
permit  interpolation  to  the  nearest  1/32  pound.   A  fluid  sample  of  at 
least  ten  pounds  was  used  as  a  basis  for  determining  the  rate  of  flow 
when  measurements  were  taken,  using  a  ten  gallon  stainless  steel  con- 
tainer to  collect  the  sample  of  fluid.   Comparison  of  the  volumetric 
and  weight  measurements  used  to  determine  flow  rate,  established  by 
weighing  a  4000  ml  fluid  sample  of  carefully  measured  density,  indicated 
the  percent  difference  to  be  0.6%.   This  difference  corresponded  to  the 
minimum  weight  discrimination  believed  achievable  on  the  Fairbanks  Morse 
Scale  (1/32  pound). 

Fluid  densities  were  measured  with  a  Westphal  balance  obtained 
from  the  W.  M.  Welch  Scientific  Company.   The  weights  were  identified  by 
serial  number  4028A  and  the  bulb  was  numbered  472  WHR.   This  balance  was 
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TABLE  2-1 
MODEL  OF  EXPECTED  FLOW  PERFORMANCE 


FLOW  TUBE  DIME  M  3  i  0  N  S  (I  M  C  H  E  S  ) 


FLOW  TURE  DIAMETER  0.0533 
TEST  SECTION  LENGTH   5.406 
ENTRANCE  LENGTH   3.997 
EXIT  LENGTH   1.33? 


FLUID  PROPEFTIES--VISCOSITY  (Cp)  1.000 

DENSITY  (GM/CC)  1.000 


EXPECTED  AMOUNT  REQUIRED 

DIFFERENTIAL  OF  FLUID  RESERVOIR 

PRESSURE  REQUIRED  PRESSURE 

(RSI  AND/  (POUNDS  AND/ 

(  F  T / S E C  )  (  G A L / M  I  N  )   IN.  OF  H G  .  )       L  I  T  E  R  S  )  (  P S  I  ) 


REYNOLDS  AVERAGE   FLOW 
NUMBER    VELOCITY  RATE 


20  0.0 
400.O 

600  ,0 

8  0  0,0 

10  0  0.0 

1200  .  0 

14  0  0.0 
16  0  0.0 

1 6  o  o .  n 

2000,0 

2  0  0  0  0.0 

4  0000.0 

6  0  0  0  0.0 
6  0  0  0  0.0 


0.5 

1.0 

1  •  5 

1.9 

2.4 

2.9 

3.4 

3.9 

4  .4 

4  .8 

48.5 

9  6.9 

145.4 

193.9 


0.0 

0-0 
CO 
0.0 

o.o 

0.0 

o.o 

0.0 

o.o 

o.o 


0.0  5 
0.10 
0.15 
0.21 

n  .26 

0.31 
0.36 
0  .41 
0  .46 
0.51 


0.11 
0  .23 
0.34 
0  .45 
0.57 
0  .68 
0  .79 
0  .91 

1.02 
1.13 

94  .30 

317.20 

6  4  4.9Q 
1.3    483.17 

1066.93 


0.3  42.71 
0.7  14  3.65 
1.0    292.05 


0  .0 


0,1 


0.1 
0.2 


0.2 


0.3 


0.3 


0  .4 


0.4 


0.5 


4.  7 


9.4 


14.1 


16.8 


0.021 


0.043 


0.064 


0.085 


0  .106 


0.128 


0.14  9 
0.170 


0.1.91. 


0.213 


2  . 1 2  7 
4.253 


6.380 


6.506 


0.1 

0.2 

0.3 

0.5 

0.6 

0.7 

0.9 

1.0 

1.1 

1.3 

111.2 

390.5 

516.2 

1379.0 


AT  AND.  ABOVE  a  REYNOLDS  NUMBER  OF   80000.0  TH=  REQUIRED  RESERVOIR 
PRESSURE  E  X  C  E  c  D  ^  THE  AVAILABLE  SUPPLY 
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initially  calibrated  by  using  distilled  water,  the  density  of  which  has 
been  established  accurately  (16).   The  balance  was  placed  on  a  smooth 
aluminum  plate  in  a  semi-enclosed  cabinet  to  minimize  environmental  dis- 
turbance.  The  calibration  was  checked  periodically  during  the  research 
and  calibration  phases  to  insure  against  calibration  shift  caused  by 
accidental  movement.   Using  this  balance,  densities  were  measured  in 
increments  of  0.0001  grams  per  cubic  centimeter  up  to  a  maximum  of  1.8999 
grams  per  cubic  centimeter. 

Temperatures  of  both  the  environment  and  the  fluids  were  taken 
using  a  Princo  15-30°C  thermometer  calibrated  in  increments  of  0.1°C  with 
a  sufficient  spread  between  marked  intervals  to  permit  interpolation 
to  the  nearest  0.05°C. 

Three  different  pressure  regulators  connecting  the  nitrogen  source 
cylinder  with  the  fluid  reservoir  were  used  depending  upon  the  range  of 
pressures  required  to  drive  the  fluid  through  the  tubes.   Regulated 
pressures  of  200-1500  pounds  per  square  inch  were  attained  using  a 
Hoke-Phoenix  0-3000  psi  regulator,  style  521  serial  number  N296C. 
Pressures  of  75-175  psi  were  obtained  through  the  use  of  a  REGO  0-400  psi 
regulator,  type  XL.   A  Victor  Equipment  Company  0-60  psi  regulator  was 
used  when  pressures  ranging  from  10-60  psi  were  required. 

Viscosity  measurements  were  taken  with  a  calibrated  Cannon- 
Fenske  capillary  tube  viscometer  number  25/J680.   The  supplied  calibra- 
tion constant  of  the  viscometer  was  0.002225  centistokes  per  second  at 
100°F,  and  0.002218  centistokes  at  210°F.   The  calibration  constant  at 
70°F  obtained  through  linear  extrapolation  of  supplied  values  was  deter- 
mined to  be  0.002227  centistokes  per  second. 

A  series  of  calibration  runs  to  check  this  value  was  conducted 
using  distilled  water  for  which  accurately  measured  data  on  viscosity 
and  density  as  a  function  of  temperature  was  readily  available  (16).   A 
computer  program  (program  number  2  in  Appendix  A)  was  then  prepared 
which  would  accept  values  of  measured  time  and  temperature  and  the  cor- 
responding known  values  of  viscosity  and  density  of  the  distilled  water. 
From  this  information  an  average  calibration  constant  was  determined  for 
the  viscometer  over  the  environmental  temperature  ranges  encountered. 
The  results  of  these  calibration. runs  are  included  in  Table  2-2.   The 
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TABLE  2-2 


RESULTS  ON  THE  CALIBRATION  OF  NO.  25/J680 
CAPILLARY  TUBE  VISCOMETER 


COMPUTED 
FLOW    FLUID   STANDARD  STANDARD   KINEMATIC  CALIBRATION 
TIME    TEMP.   DENSITY   VISCOSITY  VISCOSITY  CONSTANT 
(SECS)  (DEG.C)  (GM/CC)      <CP)       (CS)     (CS/SEC) 


431.4 
431.6 
432.8 
4  31.0 
426.0 
424.0 


21  .6 
21.6 
21.8- 
22.3 

2?. 6 
2?. 9 


0  ,9979 
0 . 9979 
0.9978 
0  .9977 
0  .9977 
0  .9976 


0.9669 
0.9669 
0.9623 
0.9510 
0 .9444 
0.9378 


0 .9690 
0 .9690 
0 . 9644 
0.9532 

0 .9466 
0 .9400 


0 .0022461 
0 .0022*51 
0 .0022283 
0 .0022116 

0 .0022220 
0 .0022170 


COMPUTED  ARITHMETIC  AVE.  CALIBRATION  CONSTANT  0.0022283 
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average  calibration  constant  of  0.002228  was  used  throughout  the  exper- 
imental phase  of  this  study  as  the  basis  for  determining  viscosity  with 
the  Cannon-Fenske  viscometer. 

The  pressure  and  differential  pressure  measuring  instruments  used 
in  these  experiments  include  a  Heise  0-2000  psi  bourdon  pressure  gage 
number  H2939R,  a  Midwest  Instrument  0-1000  psi  differential  bourdon 
pressure  gage  and  a  50  inch  glass  manometer  tube  suitable  for  safe  use 
at  pressures  corresponding  to  approximately  two  atmospheres.   Both 
bourdon  gages  are  marked  to  indicate  increments  of  5  psi  over  their 
entire  range.   The  spread  of  the  marked  increments  on  the  Heise  gage  is 
sufficient  to  permit  reliable  interpolation  to  within  1  psi  while  the 
design  of  the  Midwest  gage  is  such  that  a  reliable  interpolation  to 
within  2  psi  was  the  best  accuracy  believed  attainable. 

The  Heise  and  Midwest  gages  were  initially  calibrated  against  a 
precision  dead  weight  pressure  measuring  instrument,  Ruska  Balance 
#11516  with  piston  #B3-212,  in  the  range  from  100-1000  psi.   While  the 
Heise  gage  proved  to  be  in  exact  agreement  (within  the  1  psi  limit 
attainable  for  this  gage)  over  the  entire  range  considered,  the  corres- 
ponding readings  taken  with  the  Midwest  gage  raised  a  serious  question 
of  repeatable  reliability.   In  order  to  establish  a  statistical  base  from 
which  to  judge  the  repeatability  and  usefulness  of  this  gage,  a  series 
of  recalibration  runs  were  conducted.   The  Heise  gage  and  Midwest  gage 
were  connected  in  parallel  to  a  source  of  bottled  nitrogen  through  the 
Hoke-Phoenix  pressure  regulator  and  the  pressure  was  systematically 
raised  and  lowered  over  a  range  from  0-700  psi.   The  Heise  gage  was 
used  as  a  secondary  standard  against  which  the  performance  of  the  Midwest 
gage  could  be  measured.   The  results  of  this  calibration  are  indicated 
in  Figure  2-2.   A  mechanical  hysteresis  of  2  psi  over  this  pressure 
range  is  readily  apparent. 

While  the  statistical  scatter  was  considerable  in  the  lower  pres- 
sure ranges ,  it  was  decided  that  the  repeatability  was  sufficiently  good 
above  the  75  psi  range  to  justify  the  continued  use  of  this  gage.  How- 
ever an  allowance  providing  for  the  observed  statistical  uncertainty  has 
been  included  in  all  calculations  derived  from  values  measured  with  this 
gage.   While  it  is  recognized  to  be  the  measuring  instrument  with  the 
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least  desirable  characteristics  of  all  the  equipment  used,  it  was  the 
only  differential  pressure  gage  available  that  was  capable  of  measuring 
the  high  differential  pressures  expected  to  be  encountered  in  the 
course  of  the  experimental  investigation.   Further,  time  limitations 
would  not  permit  the  further  delays  which  would  be  required  to  obtain, 
install  and  calibrate  a  different  gage. 

CONSTRUCTION  AND  CALIBRATION  OF  THE  VARIABLE  PRESSURE  FLOW  SYSTEM 

A  specially  adapted  high  pressure  nitrogen  cylinder  was  employed 
as  the  fluid  reservoir  of  the  flow  system.   A  safe  working  pressure  of 
2200  pounds  per  square  inch  was  determined  by  tests  conducted  in  July 
of  1962.   The  capacity  of  the  cylinder  was  determined  to  be  1.78  cubic 
feet  (111.5  pounds  of  tap  water  measured  at  17.9°C).   The  exit  to  the 
nitrogen  cylinder  was  fitted  with  a  specially  designed  cap  capable  of 
receiving  nitrogen  gas  under  high  pressure  and  discharging  the  sample 
fluid  to  a  flow  tube.   A  detailed  drawing  of  the  cap  is  shown  in 
Figure  2-3. 

One  of  the  unique  features  of  the  flow  system  is  the  use  of 
specially  constructed  tubes  (detailed  drawings  of  these  tubes  are  shown 
in  Figure  2-4).   Hypodermic  tubing  with  inside  diameters  varying  from 
0.054  up  to  0.106  inches  was  obtained  from  the  C.  A.  Roberts  Company. 
Five  tubes  were  constructed  each  with  test  sections  of  approximately  100 
diameters,  entrance  lengths  of  75  diameters  and  exit  lengths  of  25 
diameters.   A  small  hole  was  drilled  at  either  end  of  each  test  section 
to  which  1/8  inch  0.  D.  stainless  steel  tubing  was  connected  through 
5/16  inch  0.  D.  copper  sleeves.   The  sleeves  were  reamed  to  snugly 
enclose  the  tubes,  thereby  providing  support  to  the  flow  tubes  as  well 
as  to  the  1/8  inch  tubing  through  which  pressure  measurements  could  be 
taken.   The  exit  end  of  the  flow  tube  discharged  to  a  larger  tube  which 
was  used  to  carry  and  deflect  the  efflux  to  a  collection  receptacle. 
The  entrance  end  of  the  tube  was  fitted  into  a  3/8  inch  Autoclave 
Engineering  plug  fitting  which  could  then  be  assembled  into  an  Auto- 
clave series  30,  90°  angle  valve.   The  entrance  rim  of  the  Autoclave 
plug  was  rounded  to  meet  the  tubes  thereby  minimizing  the  disturbance 
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25 
caused  by  an  abrupt  decrease  in  flow  diameter  at  that  particular  point 
in  the  flow  system.   Hypodermic  tubes  with  an  outside  diameter  very 
nearly  equal  to  the  inside  diameter  of  the  flow  tubes  were  used  to 
deburr  and  polish  the  flow  tubes  after  they  had  been  constructed  and 
during  their  experimental  use. 

The  length  of  the  test  section  of  the  flow  tube  was  measured  to 
better  than  1/64  inch  using  a  precision  machinist's  ruler.   However, 
since  the  internal  diameter  of  the  hydraulic  tubing  supplied  by  the  manu- 
facturer was  guaranteed  to  be  accurate  within  limits  of  approximately 
±0.002  inches,  it  was  necessary  to  determine  this  measurement  more  pre- 
cisely by  conducting  flow  runs  in  the  laminar  flow  regime.   For  this 
purpose  tap  water  was  loaded  into  the  fluid  reservoir  and  forced  through 
the  flow  tubes  at  Reynolds  numbers  no  greater  than  2,000.   A  manometer 
containing  colored  carbon  tetrachloride  which  had  a  measured  density  of 
1.557  grams  per  cubic  centimeters  was  connected  across  the  pressure 
ports  of  the  flow  tube.   Since  the  difference  in  density  of  the  mano- 
metric  fluid  and  the  tap  water  was  approximately  0.56  grams  per  cubic 
centimeters,  this  device  provided  a  very  sensitive  and  accurate  means 
of  determining  the  pressure  differential  across  the  test  section. 

A  computer  program  was  then  prepared  to  compute  the  internal 
diameter  of  the  flow  tubes  from  the  following  measured  values:   (a)  test 
section  length  in  inches;   (b)  fluid  density  in  grams  per  cubic  centi- 
meter; (c)  fluid  viscosity  as  measured  by  the  Cannon-Fenske  viscometer 
in  seconds;  (d)  volume  of  fluid  flow  in  milliliters;   (e)  the  time  of 
fluid  flow  in  seconds;  and   (f)  pressure  in  inches  of  manometric  fluid. 
The  basis  for  the  computation  is  the  well  known  Hagen-Poiseuille  law: 

Q  =  tt(P   -  Pn)  g  R4/8yL  [2-3] 

o    1   c 

which  appears  in  computer  program  #3  in  Appendix  A  in  the  form  of: 

D  =  [(128yQL)/(TrAPgc)r  [2-4] 

The  results  of  these  calibration  runs  are  summarized  in  Table  2-3.  Com- 
puter listings  for  each  of  the  five  flow  tubes  calibrated  in  this  manner 
appear  in  Appendix  B. 
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TABLE  2-3 
SUMMARY  OF  FLOW  TUBE  DIMENSIONS 


Manufacturers 

Test 

Flow  Tube 

Supplied 

Calibrated 

Section 

Entrance 

Exit 

Gage 

Diameter 

Diameter 

Length 

Length 

Length 

Number 

(Inches) 

(Inches) 

( Inches) 

( Inches) 

( Inches ) 

15 

0.054 

0.05331 

5.406 

3.99 

1.33 

13 

0.071 

0.07158 

7.109 

5.325 

1.78 

12 

0.085 

0.08548 

8.500 

6.375 

2.13 

11 

0.094 

0.09264 

9.460 

7.05 

2.35 

10 

0.106 

0.10590 

10.609 

7.95 

2.65 

In  addition  to  providing  a  means  by  which  the  internal  diameter 
of  the  flow  tubes  could  be  accurately  measured,  the  results  of  these 
tests  indicated  that  data  of  reasonable  consistency  could  be  obtained 
in  the  laminar  flow  regime  using  the  measurement  techniques  and  instru- 
ments employed. 

Performance  characteristics  of  the  tubes  was  next  tested  in  the 

turbulent  flow  regime  corresponding  to  a  range  of  Reynolds  numbers  from 

4      5 

10  to  10  .   While  it  would  have  been  desirable  to  measure  these  char- 
acteristics for  each  of  the  five  tubes,  time  limitations  permitted 
detailed  examination  of  only  two  tubes,  number  15  and  13.   These  calibra- 
tion experiments  were  conducted  over  the  same  turbulent  flow  regime 
using  two  fluids,  military  jet  fuel  (JP4)  and  tap  water.   The  flow 
behavior  of  each  of  these  fluids  was  thereby  independently  reproduced. 

There  were  three  major  objectives  to  be  accomplished  during 
these  performance  tests: 

(a)  To  determine  the  interior  roughness  characteristics  of 
the  'flow  tubes . 

(b)  To  determine  the  effects,  if  any,  of  the  pressure  ports 
placed  at  either  end  of  the  test  section  over  the  range  of  Reynolds 
stresses  to  which  the  fluids  would  be  subjected  during  research  experi- 
ments. 
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(c)   To  observe  the  performance  of  the  Midwest  differential 
pressure  gage  as  a  measuring  device  under1  dynamic  turbulent  flow  condi- 
tions . 

It  might  be  well  to  recall,  in  connection  with  this  latter 
objective,  that  the  characteristics  of  the  Midwest  differential  pres- 
sure gage  were  originally  determined  by  a  calibration  technique  which 
produced  a  static   pressure  differential  across  the  two  inlet  ports  of 
the  gage.   Further,  there  was  a  clearly  defined  mechanical  hysteresis 
associated  with  this  gage  which  appeared  when  static  pressures  were 
being  systematically  increased  and  decreased.   This  observed  static 
behavior  raised'  serious  questions  of  the  gage's  ability  to  yield  reli- 
able pressure  measurements  under  conditions  of  turbulent  fluid  flow 
which  is  characterized  by  violent  pressure  fluctuations. 

It  was  recognized,  of  course,  that  the  effects  of  the  pressure 
ports  could  only  be  judged  if  the  interior  roughness  characteristics 
failed  to  conform  with  well  established  flow  behavior  data.   The  dilemma 
just  described  is  not  uncommon  to  scientific  research:   Experimental 
results  cannot  be  achieved  without  measurement — yet  to  measure  is  to 
disburb.   Further,  whether  a  meaningful  discrimination  between  the 
first  two  objectives  could  be  achieved  would  largely  depend  upon  the 
performance  of  the  differential  pressure  gage. 

There  were,  however,  some  criteria  by  which  the  results  of  these 
tests  could  be  judged: 

(a)  The  results  of  the  calibration  tests  on  the  differential 
pressure  gage  suggested  a  decrease  in  reliability  at  the  lower  pressure 
ranges.   Further,  at  pressure  readings  of  less  than  100  psi,  the  number 
of  significant  digits  available  was  reduced  by  one. 

(b)  Any  abnormal  effects  on  flow  behavior  resulting  from  the 
pressure  ports  could  be  expected  to  become  more  pronounced  as  the  flow 
rates  were  increased. 

(c)  The  expected  minimum  friction  factor  at  any  particular 
flow  rate  would  correspond  to  the  friction  factor  established  for 
"hydraulically  smooth"  pipes.   The  expected  maximum  would  not  be  likely 
to  exceed  the  friction  factor  relationship  established  for  drawn  tubing. 
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For  the  15  gage  flow  tube  with  a  0.053  inch  I.D.,  the  maximum  expected 
relative  roughness  was  approximately  0.002.   The  corresponding  value 
for  the  13  gage  tube  was  approximately  0.001. 

Using  these  criteria,  the  results  of  the  turbulent  flow  tests 
as  illustrated   in  Figure  2-5  suggest  the  following  conclusions: 

(a)  The  reliability  of  pressure  measurements  taken  with  the 
Midwest  gage  is  substantially  reduced  at  differential  pressure  readings 
below  75  psi. 

(b)  The  relative  roughness  of  the  flow  tubes,  including  any 
contributing  effects  of  the  pressure  ports,  is  between  1/5  and  1/10  the 
relative  roughness  associated  with  drawn  tubing. 

(c)  No  abnormal  flow  behavior  can  be  specifically  identified 
with  the  pressure  ports  over  the  range  of  shear  rates  examined  during 
the  test. 

Hence,  it  was  concluded  that  the  equipment  and  flow  system  design 
were  entirely  satisfactory  for  experimental  studies  on  drag  reduction. 
A  computer  write  out  of  information  used  to  construct  Figure  2-5  is 
contained  in  Appendix  C.   Computer  program  number  4 ,  used  to  produce 
this  information  is  contained  in  Appendix  A. 

OPERATION  OF  FLOW  SYSTEM 

Once  a  sample  had  been  loaded  into  the  fluid  reservoir  a  choice 
of  the  flow  tube  to  be  used  and  the  range  of  the  flow  regime  to  be 
investigated  was  required  in  order  to  determine  appropriate  measurement 
equipment.   Generally,  the  Midwest  differential  pressure  gage  was 
employed  to  measure  the  pressure  drop  across  the  test  section  at  high 
Reynolds  numbers  and  the  flow  rate  was  measured  on  the  Fairbanks  Morse 
Scale  by  timing  a  predetermined  weight   of  fluid.   The  mercury  manometer 
was  employed  to  measure  the  test  section  pressure  differential  at  the 
lower  Reynolds  numbers  and  the  flow  of  a  predetermined  volume   of  fluid 
(up  to  a  maximum  of  4000  milliliters)  was  timed  to  measure  flow  rate. 

The  computer  flow  model  commented  upon  earlier  in  this  Chapter 
was  used  as  a  guide  for  these  determinations.   By  referring  to  Table  2-1, 
for  example,  if  the  flow  of  pure  tap  water  through  the  15  gage  tube  at  a 
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Reynolds  number  in  the  vicinity  of  60,000  was  to  be  observed,  the  pres- 
sure regulator  from  the  nitrogen  cylinder  was  set  at  800  pounds  per 
square  inch  and  the  program  indicated  that  the  flow  of  at  least  15 
pounds  (rounded  to  the  nearest  5  pounds)  of  the  fluid  should  be  timed 
in  order  to  insure  that  the  timed  run  would  consume  at  least  100  sec- 
onds.  Therefore,  a  series  of  runs  could  be  preplanned  to  maximize  the 
efficient  use  of  fluid  available  for  investigation. 

Once  the  flow  regime  to  be  investigated  had  been  determined  and 
the  appropriate  measurement  equipment  assembled,  it  was  next  necessary 
to  purge  the  tubing  leading  from  the  pressure  ports  of  the  flow  tubes 
to  the  pressure  gages.   This  was  accomplished  by  loosening  the  fittings 
at  the  entrance  to  the  gages ,  plugging  the  exit  end  of  the  flow  tube 
and  permitting  the  fluid  to  flow  through  the  pressure  gage  line  until 
a  free  flow,  clear  of  gas  bubbles  was  observed.   Next,  a  500  milliliter 
sample  of  the  fluid  was  withdrawn  from  the  reservoir  from  which  the 
basic  properties  of  temperature,  density  and  viscosity  were  determined. 
An  average  measurement  of  viscosity  was  obtained  on  the  sample  by 
timing  three  to  five  consecutive  trial  runs  in  the  Cannon-Fenske  capil- 
lary tube  viscometer.   Since  there  was  not  a  means  provided  to  maintain 
the  reservoir  fluid  at  a  constant  temperature,  the  fluid  properties  of 
the  reservoir  mixture  were  measured  every  20  to  30  minutes  during  the 
course  of  experimentation.   Variations  in  these  properties  for  a  given 
reservoir  load  were  generally  found  to  be  negligible  since  the  maximum 
time  required  to  empty  the  cylinder  rarely  exceeded  45  minutes  once 
experiments  were  begun.   The  cylinder  used  with  the  Westphal  balance  and 
the  capillary  tube  viscometer  were  each  cleaned  and  dried  immediately 
after  use. 

Prior  to  a  test  run,  the  weights  on  the  Fairbanks  Morse  Scale 
were  adjusted  to  permit  two  to  four  pounds  of  the  test  fluid  to  flow 
before  the  stop  watch  was  started  and  measurements  were  begun.   This  pro- 
cedure allowed  a  sufficient  amount  of  time  for  the  Autoclave  90°  valve 
to  be  fully  opened  and  permitted  the  fluid  flow  pattern  to  develop  and 
stabilize  before  measurements  were  taken.   When  the  arm  of  the  Fairbanks 
Morse  Scale  returned  to  equilibrium,  the  stop  watch  was  started  and  the 
predetermined  weights  were  added  to  the  balance  arm.   Measurements  of 
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gage  pressure  at  the  upstream  end  of  the  test  section  and  differential 
pressure  across  the  test  section  as  well  as  the  weight  of  fluid  to  be 
timed  were  recorded  on  a  data  sheet.  When  enough  fluid  entered  the 
collection  container  on  the  scale  to  cause  the  arm  to  return  to  equil- 
ibrium, the  stop  watch  was  stopped  and  the  flow  valve  was  closed.  The 
temperature  of  the  fluid  was  then  taken  and  recorded  together  with  the 
elapsed  time  for  the  run. 

ANALYSIS  OF  FLOW  DATA 

One  basis  for  determining  the  effectiveness  of  a  polymer  would 
be  to  relate  the  friction  factor  observed  for  a  fluid  with   such  an 
additive  to  the  friction  factor  observed  for  a  purely  viscous  fluid 
flowing  at  the  same  Reynolds  number.   The  factors  affecting  the  Reynolds 
number- - ie. ,  density,  average  fluid  velocity,  tube  diameter  and  fluid 
viscosity--could  all  be  accurately  determined  to  three  or  more  signifi- 
cant figures.   However,  the  friction  factor  varies  as  a  direct  function 
of  pressure  which,  when  measured  with  the  Midwest  differential  pressure 
gage,  could  be  determined  to  two  significant  figures  at  best.   Informa- 
tion indicating  the  scatter  in  pressure  readings  expected  over  the  wide 
range  of  pressures  investigated  was  available  from  the  calibration  runs 
performed  on  the  gage  (see  Figure  2-2).   This  scatter  includes  an  allow- 
ance for  both  the  observed  mechanical  hysteresis  as  well  as  the  random 
scatter  associated  with  other  limitations  of  the  instrument  and  observer, 
This  information  was  used  to  compute  error  bounds  for  each  value  of 
friction  factor,  assuming  that  the  measurements  taken  with  the  differen- 
tial pressure  gage  represented  the  least  reliable  data  and  would 
contribute  the  most  significant  error. 

At  the  conclusion  of  a  series  of  runs,  data  cards  were  prepared 
using  the  information  taken  during  each  run.   These  were  used  with 
Computer  Program  number  4  in  Appendix  A  to  analyze  the  test  data. 

PREPARATION  OF  POLYMER  SOLUTIONS 


The  pioneering  research  undertaken  by  Fabula  and  Hoyt  (6)  at  the 
Naval  Ordnance  Test  Station  during  the  late  1950' s  and  early  1960's 
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was  instrumental  in  identifying  a  relatively  large  number  of  water 
soluble  polymers  capable  of  inducing  drag  reduction  effects  in  extremely 
dilute  aqueous  solutions.   However,  an  equally  extensive  program 
directed  toward  the  identification  of  additives  capable  of  producing 
drag  reduction  effects  in  dilute,  non-polar  hydrocarbon  solutions  has 
never  been  undertaken.   In  fact,  prior  to  Shin's  investigation  (23), 
Toms  original  research  (25)  appears  to  be  unique,  not  only  because  it 
represents  the  first  written  account  of  drag  reduction  associated  with 
dilute  Newtonian  solutions,  but  also  because  his  experiments  involved  a 
non-aquecus  solvent.   Shaver  and  Merrill  worked  with  concentrated   solu- 
tions of  polyisobutylene  in  cyclohexane  (22).   During  World  War  II,  an 
apparent  reduction  of  flow  resistance  was  noted  when  gasoline  was 
thickened   with  napalm  for  use  in  flame  throwers  (1).   However,  in 
contrast  to  the  dilute  solutions  studied  by  Toms  and  Fabula.  the  final 
solutions  involved  in  these  isolated  examples  had  markedly  different 
fluid  characteristics  than  the  original  solvents.   The  fact  that  poly- 
isobutylene had  been  shown  to  be  soluble  in  a  pure  non-polar  hydro- 
carbon solvent  and  the  fact  that  it  apparently  displayed  anomalous 
behavior  in  concentrated  solutions  clearly  suggested  that  further  study 
of  its  behavior  in  dilute  solutions  might  be  productive.   Additionally, 
it  seemed  appropriate  to  test  the  more  common  water  soluble  polymers  for 
solubility  in  jet  fuel.   For  this  reason,  samples  of  several  polymers 
known  to  exhibit  viscoelastic  properties  were  initially  requested  and 
obtained  from  various  industrial  sources.   While  most  of  these  polymers 
were  known  to  be  soluble  in  water  and  sea  water,  the  solubility  of  each 
was  tested  in  the  jet  fuel  for  possible  further  examination  if  time  would 
permit.   A  summary  of  the  polymer  samples  and  the  sources  from  which  they 
were  obtained  is  contained  in  Table  2-U. 

Several  50  weight  part  per  million  solutions  in  jet  fuel  were  pre- 
pared, each  containing  one  of  these  polymers ,  and  all  were  observed  for 
several  days.   Each  solution  was  prepared  by  carefully  weighing  a  small 
sample  of  the  polymer  (between  34  and  38  milligrams)  on  a  precalibrated 
precision  Metier  analytic  balance  and  then  adding  the  volume  of  jet  fuel 
of  carefully  measured  density  that  would,  produce  a  concentration  of  50 
weight  parts  per  million.   These  samples  were  placed  in  a  one  liter 


TABLE  2-4 
LIST  OF  POLYMERS  OBTAINED  AND  INDUSTRIAL  SOURCE 
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Company 


Product 


Archer  Daniels  Midland 
Minneapolis,  Minn, 

En jay  Company 
New  York,  N.Y. 

Union  Carbide  Corporation 
New  York,  N.Y. 


Hercules  Powder  Company 
Wilmington,  Del. 

Stein,  Hall  S  Co. ,  Inc. 
New  York,  N.Y. 


Polysac.chori.de 
ADM  7079 

Polyisobutylene 

L80,  L140  and  L200 

Polyox  WSR  P301, 
WSR  P205, 
N3000  and  N750 

Reten  210,  Al  and  A5 


Jaguar  J2S1 
Polyhal  27 


bottle  and  were  periodically  tested  for  solubility  by  visual  observation. 
None  of  the  water  soluble  polymers  appeared  to  be  soluble  in  jet  fuel  when 
mixed  by  natural  diffusion  techniques  at  room  temperature.   Although  mix- 
ing by  diffusion  can  be  a  very  slow  and  difficult  process  for  many  poly- 
mers, it  did  not  appear  that  more  vigorous  artificial  methods  should  be 
employed  in  view  of  the  well  known  susceptibility  of  most  polymers  to 
mechanical,  thermal  and/or  chemical  degradation.   On  the  other  hand,  the 
three  samples  of  polyisobutylene  prepared  in  one  liter  bottles  at  con- 
centrations of  50  weight  parts  per  million  appeared  to  be  completely  dis- 
solved through  natural  diffusion  when  left  standing  at  room  temperature 
for  approximately  36  hours.   The  preparation  of  more  concentrated  samples 
was  found  to  be  somewhat  slower.   However,  even  a  6%  weight  solution  of 
the  highest  molecular  weight  polymer  prepared  (VISTANEX  L-200)  was 
completely  dissolved  as  judged  by  visual  inspection  of  such  factors  as 
homogeneity,  clarity,  etc.,  in  less  than  two  weeks.   The  speed  with  which 
the  mixing  occurred  appears  to  be  more  a  function  of  the  total  quantity 
of  polymer  present  in  a  sample  than  the  weight  concentration  of  the 
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sample  being  prepared.   That  is,  a  3%  solution  prepared  in  500  ml  size 
samples  dissolved  more  rapidly  than  a  3%  solution  prepared  in  1000  ml 
size  batches. 

The  concentrated  polymer  solutions  were  developed  to  facilitate 
the  preparation  of  the  50  weight  part  per  million  solutions  to  be 
loaded  into  the  fluid  reservoir  for  flow  experiments.   It  was  found  that 
these  solutions  once  the  polymer  was  fully  dissolved,  mixed  very  readily 
with  the  pure  solvent  to  form  the  desired  dilute  concentration.   Thus,  by 
pouring  a  solution  containing  2.267  grams  of  the  polymer  into  a  20  gallon 
stainless  steel  tank  and  then  adding  enough  jet  fuel  to  bring  the  total 
weight  to  100  pounds,  a  50  wppm  sample  could  be  prepared  in  a  very  few 
minutes.   In  an  effort  to  insure  that  all  of  the  polymer  originally 
prepared  in  the  one  liter  bottle  was  actually  included  in  the  100  pound 
solution,  the  bottle  was  flushed  several  times  with  the  pure  jet  fuel 
which  was  used  to  form  the  final  solution. 

Since  the  preparation  of  large  quantities  of  even  dilute  concen- 
trations can  become  a  problem,  the  procedure  adopted  in  this  study  is 
certainly  recommended.   Further,  certain  physical  properties  which  are 
impossible  to  describe  quantitatively  such  as  homogeneity,  purity,  etc. 
become  more  difficult  to  control  as  the  size  of  the  sample  increases. 

Once  the  desired  dilute  solution  was  formed,  it  was  gently  poured 
into  the  fluid  reservoir  at  a  rate  of  about  1/2  gallon  per  minute  to 
minimize  any  effects  of  mechanical  degradation  or  depolymerization  prior 
to  actual  flow  through  the  test  section. 


CHAPTER  III 
EXPERIMENTAL  RESULTS  AND  DISCUSSION 

All  drag  reduction  experiments  were  accomplished  by  observing  the 
flow  behavior  of  solutions  containing  50  weight  parts  per  million  of  the 
high  molecular  weight  Vistanex  L-200  in  military  jet  fuel.   The  two  flow 
tubes  with  the  smallest  diameters  were  used:  number  13  with  an  0.0716 
inch  I.  D.  and  number  15  with  an  0.0533  inch  I.  D.   The  flow  system  was 
primarily  operated  in  the  turbulent  regime  corresponding  to  a  range  of 
Reynolds  numbers  from  approximately  10,000  to  100,000,  although  a  few 
observations  were  also  made  on  laminar  flow.   Values  of  observed  friction 
factor  and  Reynolds  number  were  calculated  for  each  run  using  the  same 
computer  program  (number  4  in  Appendix  A)  as  was  used  to  analyze  flow 
data  pertaining  to  the  pure  solvents. 

The  results  of  the  turbulent  flow  experiments  are  illustrated  in 
Figure  3-1.   (The  data  is  presented  in  Appendix  D.)   The  friction  factor 
versus  Reynolds  number  relationship  found  to  exist  for  purely  viscous 
fluids  flowing  in  the  turbulent  regime,  as  well  as  an  extension  of 
laminar  flow  relationships  have  been  included  in  order  to  facilitate 
visual  comparison.   Error  bands  have  been  included  on  all  data  calculated 
from  pressure  measurements  taken  with  the  Midwest  D.P.I,  gage  in  order  to 
illustrate  the  statistical  uncertainty  observed  at  the  time  the  gage  was 
calibrated.   Flow  lines  have  been  drawn  which  place  greatest  reliability 
on  the  data  points  calculated  from  pressure  measurements  taken  with  the 
mercury  manometer.   These  data  points  do  not  have  error  bands. 

In  examining  Figure  3-1  it  should  be  noted  that: 

(a)  At  any  given  Reynolds  number,  the  drag  reduction  effect 
appears  to  be  more  pronounced  as  the  internal  diameter  of  the  flow  tube 
increases . 

(b)  Over  the  range  of  Reynolds  stresses  covered  by  these 
experiments ,  the  drag  reduction  effect  appears  to  first  increase  then 
decrease  with  increasing  Reynolds  numbers. 

The  former  occurrence  is  well  known  and  has  been  reported  since 
Toms  original  findings  as  the  "diameter  effect."  Its  importance,  how- 
ever, is  often  overlooked.   It  is  because  of  this  effect  that  dynamic 
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Fig.  3-1.   Log-Log  Plot  of  Observed  Friction  Factor  vs.  Reynolds 
Number  for  the  Flow  of  50  wppm  PIB  in  Jet  Fuel 
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similarity  in  the  turbulent  flow  regime  is. lost.   The  Reynolds  Number 
Similarity  Law  for  turbulent  flow  of  fluids  no  longer  applies  to  this 
unique  class  of  fluids.   On  the  other  hand,  this  anomalous  behavior 
does  not  occur  in  the  laminar  region  of  fluid  flow.   A  friction  factor 
equal  to  or  slightly  higher  than  the  relationship  which  correlates  the 
Fanning  friction  factor  with  the  Reynolds  number  in  the  laminar  regime : 

f  =  16/NRe  [3-1] 

was  found  to  apply  to  this  c-lass  of  fluids.   (This  data  is  not  illus- 
trated on  Figure  3-1  which  has  been  scaled  so  as  to  permit  graphical 

...  .  .  4 

amplification  of  the  more  important  region  from  N   >10  .   However,  the 

statement  can  be  verified  by  reference  to  the  complete  data  included  in 
Appendix  D. ) 

The  latter  occurrence  is  somewhat  more  unique.   In  general,  it 
has  been  found  that  the  drag  reduction  effect  becomes  more  pronounced 
at  increasing  flow  rates.   Many  experimenters  (6),  (18),  (19),  (23), 
have  observed  that  there  is  a  critical  flow  rate  beyond  which  the  drag 
reducing  effectiveness  of  an  additive  appears  to  deteriorate.   Precisely 
what  determines  this  critical  flow  rate  has  not  yet  been  determined. 
The  critical  flow  rate  appears  to  be  very  dependent  upon  polymer  con- 
centration, temperature,  the  particular  polymer  used  etc.,  although 
these  observations  are  not  confirmed  in  this  study  since  only  one  poly- 
mer at  one  concentration  was  examined  at  ambient  temperatures.   This 
decrease  in  friction  reducing  effectiveness  was  found  to  occur  at 
around  a  Reynolds  number  of  50,000.   It  is  believed  to  be  indicative 
of  severe  mechanical  degradation. 

At  this  point  in  the  experimental  program,  it  was  believed  that 
greater  insight  into  the  drag  reduction  phenomena  could  be  achieved  by  a 
study  of  this  apparent  degradation  than  would  be  obtained  by  continuing 
on  the  plan  originally  established  for  the  study.   This  departure 
appeared  appropriate  for  a  number  of  reasons : 

(a)   The  recent  work  of  Shin,  (23),  Hershey,  (8),  Rodriguez 
et  al.  (20)  had  served  to  confirm  the  drag  reducing  behavior  of  dilute 
concentrations  of  polymers  in  hydrocarbon  solvents.   The  effects  of 
varying  concentrations  and  flow  tube  diameters  were  well  documented. 


38 

(b)   Any  practical  application  of  the  phenomenon  in  the 
field  of  military  logistics  as  envisioned,  must  consider  the  limita- 
tions imposed  by  mechanical  degradation. 

This  is  not  to  suggest  that  these  limitations  cannot  be  advanta- 
geously employed.  Since  the  particular  additive  employed  in  this  study 
consists  of  a  hydrocarbon  polymer  of  paraffinic  structure  (26),  some 
form  of  artificial  degradation  after  fuel  transmission  has  been 
accomplished  might  well  be  preferred  as  a  means  of  overcoming  fuel 
quality  deterioration. 

Two  methods  of  artificially  imposing  mechanical  degradation  were 
employed:   (a)   The  first  method  involved  the  reuse  of  a  solution  which 
had  already  been  subjected  to  experimentation.   The  degradation  imposed 
in  this  manner  might  well  be  representative  of  the  effects  to  be  expected 
as  a  result  of  long  distance  transmission;   (b)   The  second  method 
involved  a  modification  of  the  flow  system.   The  modification  intro- 
duced a  means  of  imposing  a  very  high  localized  shear  point  prior  to 
the  flow  tube.   It  will  be  recalled  that  care  was  exercised  in  the 
original  design  of  the  system  to  eliminate  the  occurrence  of  excessive 
fluid  shearing  prior  to  the  point  at  which  the  fluid  entered  the  flow 
tube.   Even  the  entrance  to  the  flow  tube  was  tapered  to  minimize  the 
pressure  loss  at  that  location.   Indeed,  when  the  15  gage  flow  tube  was 
used  to  observe  fluid  flowing  at  a  Reynolds  number  of  100,000  in  the 
tube,  the  maximum  Reynolds  number  achieved  at  any  other  location  in  the 
system  was  approximately  15,000. 

The  modification  for  the  second  method  of  imposing  mechanical 
degradation  was  accomplished  by  inserting  a  needle  valve  into  the  system 
just  upstream  of  the  flow  tube.   A  very  high  pressure  drop  could  be 
imposed  across  this  needle  valve  to  produce  an  extremely  high  localized 
Reynolds  number.   The  degradation  imposed  in  this  fashion  would  tend  to 
simulate  the  mechanical  devices  often  found  in  a  complete  closed  loop 
flow  system  which  produce  high  localized  Reynolds  stresses  such  as 
centrifugal  pumps,  etc.   This  system  is  referred  to  as  the  "modified  flow 
system." 
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The  results  of  these  experiments  are  illustrated  in  Figure  3-2. 
The  mechanical  degradation  imposed  by  both  methods  employed  is  clearly 
apparent.   The  more  severe  form  would  appear  to  accompany  the  use  of 
the  needle  valve  in  the  flow  system.   The  polymer  appears  to  be  so 
degraded  that  the  fluid  no  longer  exhibits  any  substantial  drag  reduc- 
ing characteristics.   Even  the  solution  used  for  the  first  time  through 
this  arrangement  (labeled  as  "first  pass  through  modified  flow  system") 
appears  to  be  substantially  more  degraded  than  the  solution  being 
recycled  through  the  flow  system  for  the  second  time.   These  results  sug- 
gest that  flow  system  design  can  have  a  very  important  effect  on  the 
performance  of  a  drag  reducing  additive.   That  is,  if  maximum  performance 
from  a  drag  reducing  additive  is  to  be  achieved,  mechanical  devices 
which  produce  high  localized  Reynolds  stresses  should  be  avoided.   Addi- 
tives should  be  injected  down  stream  from  such  a  device. 

The  mechanical  degradation  associated  with  either  method  appears 
to  have  permanently  destroyed  some  of  the  drag  reducing  potential  of  the 
additive.   That  is,  the  friction  factors  observed  for  the  solution  which 
was  passed  through  the  system  for  the  fourth  time  were  consistently 
higher  than  those  observed  during  even  the  third  pass. 

Further,  there  would  appear  to  be  a  striking  similarity  between 
the  effects  on  friction  drag  resulting  from  the  use  of  a  degraded  solu- 
tion and  what  has  been  termed  "the  diameter  effect."   This  suggests  that 
the  diameter  effect  is  in  fact  more  properly  identified  as  a  degradation 
effect.   That  is,  the  friction  reducing  capability  of  a  solution  is  not 
necessarily  decreased  by  virtue  of  fluid  flow  performed  in  smaller  tubes. 
Rather,  the  mechanical  degradation  has  been  increased.   This  is  not 
unexpected  since  the  mechanical  shear  applied  to  a  fluid  flowing  at  a 
given  Reynolds  number  is  substantially  increased  as  tube  diameter  is 
decreased.   Nonetheless,  the  friction  reducing  characteristics  of  the 
additive,  even  in  a  degraded  state,  continue  to  be  pronounced  when  com- 
pared to  the  flow  behavior  of  the  pure  solvent. 

Finally,  note  the  trend  of  the  flow  behavior  patterns  developing 
at  Reynolds  numbers  between  10,000  and  50,000  as  illustrated  in  Fig- 
ure 3-2.   The  friction  factor  curve  produced  from  experiments  using  a 
fresh  50  wppm  polyisobutylene  solution  appears  to  approach  an  extension 
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of  the  laminar  flow  curve.   On  the  other  hand,  the  most  highly  degraded 
solutions  almost  coincide  with  the  turbulent  friction  factor  curve. 

These  characteristics  were  also  observed  in  a  qualitative  sense 
during  the  flow  experiments.   During  the  flow  of  tap  water,  pure  jet 
fuel  or  jet  fuel  with  highly  degraded  polyisobutylene  at  well  developed 
turbulent  flow  rates,  a  very  marked  pressure  fluctuation  was  observed 
which  was  indicated  by  a  pronounced  "flutter"  of  the  Heise  gage  pointer 
and  manometric  fluid.   This  flutter  could  be  interpreted  as  variations 
in  the  absolute  pressure  as  measured  at  the  wall  of  the  tube  correspond- 
ing to  the  turbulent  fluctuations  known  to  exist.   On  the  other  hand, 
as  the  less  degraded  solutions  were  used  there  was  a  noticeable  decrease 
in  the  magnitude  of  the  flutter  associated  with  each  pressure  indicator. 
When  the  freshly  prepared  solution  was  used,  no  noticeable  flutter 
could  be  visually  detected.   It  is  recognized  that  the  measuring  instru- 
ments employed  did  not  possess  the  sensitivity  required  to  reduce  these 
qualitative  observations  to  quantitative  data.   Nonetheless,  the 
general  pattern  just  described  could  be  visually  detected  even  if  accu- 
rate measurements  could  not  be  taken  and  the  qualitative  significance  of 
these  observations  should  not  be  ignored. 

A  practical  explanation  of  this  observed  behavior  might  be  found 
in  different  modes  of  flow  of  the  solvent  and  solute.   Although  the  flow 
of  the  solvent  is  known  to  be  turbulent  at  the  Reynolds  numbers  investi- 
gated, little  is  known  concerning  the  mode  of  flow  of  the  dissolved 
solute.   It  would  appear  that  the  combined  effect  of  solvent  and  solute 
flowing  at  a  given  Reynolds  number  produces  a  unique  flow  regime,  the 
nature  of  which  depends  upon  the  relative  (or  superimposed)  contribution 
of  each.   The  flow  characteristics  of  the  fresher,  less  degraded  solu- 
tion, least  resemble  the  turbulent  characteristics  of  the  pure  solvent. 
The  inference  is  that  the  pure  solute  contributes  a  stabilizing  influence 
more  typical  of  laminar  flow,  producing  a  "pseudo-laminar"  regime  when 
the  effect  of  the  solute  is  found  to  dominate.   The  maximum  stabilizing 
contribution  of  the  solute  could  be  expected  to  occur  at  what  has  been 
called  the  "optimum"  concentration.   At  concentrations  less  than 
optimum,  the  turbulent  contribution  of  the  solvent  begins  to  dominate. 
At  concentrations  greater  than  optimum  the  contribution  of  the  solute 
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continues  to  dominate.   However,  because  of  the  thickening  action  of 
the  solute,  the  viscous  drag  could  be  expected  to  be  somewhat  greater 
than  that  obtained  at  an  optimum  concentration. 

Mechanical  degradation  occurring  during  a  flow  experiment  can 
be  expected  to  alter  the  relative  influence  of  the  solute  by  effectively 
changing  the  original  concentration.   Thus,  a  highly  degraded  solution 
originally  containing  50  wppm  of  solute  could  be  expected  to  behave  like 
a  much  weaker  undegraded  solution.   Although  this  explanation  does  not 
attempt  to  describe  the  mechanism  responsible  for  drag  reduction 
associated  with  the  solute ,  it  suggests  that  any  attempt  to  correlate 
properties  of  solutions  which  exhibit  drag  reduction  characteristics 
must  consider  the  very  important  influence  of  mechanical  degradation. 

Unfortunately,  until  such  a  correlation  is  found,  there  is  no 
mathematical  device  through  which  the  performance  of  a  given  flow  system 
can  be  predetermined  since  the  dynamic  similarity  laws  of  purely 
viscous  fluids  do  not  appear  to  apply  to  this  unique  class  of  fluids. 
However,  it  can  be  shown  that  the  rate  of  fluid  flow  in  a  given  flow 
system  is  proportional  to  the  reciprocal  of  the  square  root  of  the 
friction  factor.   That  is, 

Q    =  Q   (f  /f  )2  [3-2] 

ve    v   v  ve 

where:   Q  represents  flow  rate; 

f  represents  friction  factor; 

v  refers  to  a  purely  viscous  fluid;  and, 

ve  refers  to  a  viscoelastic  fluid. 

Although  the  precise  value  to  be  assigned  to  f   is  not  known, 

ve 

this  relationship  is  useful  in  assessing  the  potential  benefits  to  be 
derived  from  the  drag  reduction  effect  during  tactical  military  replen- 
ishment.  For  example,  a  65%  reduction  in  friction  factor  from  a  value 
of  0.00517  to  0.00181  was  the  highest  achieved  in  this  study.   The  capac- 
ity of  the  flow  system  to  transmit  fluids  was  thereby  increased  by 
approximately  69%.   Reductions  in  friction  factor  of  80%  are  not 


uncommon  (6,  23).   If  this  reduction  could  be  achieved  at  a  Reynolds 
number  of  10  ,  the  capacity  of  the  system  would  be  increased  by  approx: 
mately  224%.   In  other  words,  the  time  required  to  transmit  fuel  under 
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these  circumstances  could  be  better  than  halved.   Even  the  solution 
passed  through  the  system  for  the  fourth  time  during  this  study 
exhibited  a  28%  decrease  in  frictional  drag  at  a  moderate  Reynolds  num- 
ber of  56,000.   The  flow  rate  was  thereby  increased  some  17.5%. 

Time  did  not  permit  the  evaluation  of  quality  deterioration 
introduced  by  the  addition  of  50  wppm  of  polyisobutylene  in  military 
jet  fuel.   It  was  always  recognized  that  the  application  of  this  con- 
cept to  military  logistics  could  not  proceed  until  this  very  important 
phase  had  been  accomplished.   However,  the  benefits  of  achieving  sub- 
stantial drag  reduction  with  relatively  small  amounts  of  additives  are 
believed  to  be  sufficiently  attractive  to  warrant  this  further  effort. 


CHAPTER  IV 
CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 

The  results  of  the  experimental  study  on  the  flow  behavior  of 
polyisobutylene  dissolved  in  jet  fuel  indicate  the  following  conclusions 
and  recommendations  for  future  work . 

1.   Performance  of  Flow  System  The  variable  pressure,  single  pass 
flow  system  designed  and  constructed  for  use  in  this  study  was  found  to  be 

particularly  well  suited  for  experimental  studies  on  viscous  drag  reduc- 

5 
tion.  Tlow  behavior  up  to  a  Reynolds  number  of  10   could  be  observed  usm| 

the  apparatus  as  constructed  for  this  study,  yet  fluid  consumption  could 

be  minimized  by  the  proper  selection  of  flow  tubes.   The  single  pass  type 

test  system  facilitated  the  study  of  shear  degradation  since  the  fluid 

makes  only  one  pass  through  the  test  section  so  that  exposure  of  the 

fluid  to  shear  stress  could  be  closely  controlled.   Also,  this  type  of 

system  is  quite  similar  to  that  used  in  tactical  military  refueling 

operations . 

The  apparatus  was  found  to  perform  well  in  both  the  laminar  and 

turbulent  flow  regimes  when  tested  with  liquids  of  known  fluid  properties. 

No  abnormal  flow  behavior  was  identified  with  the  pressure  ports  of  the 

smaller  tubes  (numbers  15  and  13)  over  the  range  of  shear  rates  examined 

during  the  calibration  tests.   The  larger  tubes  were  not  used  for  drag 

reduction  flow  tests.   However,  from  an  examination  of  the  data  taken  to 

determine  the  internal  diameter  of  the  larger  tubes,  it  would  appear 

that  the  pressure  ports  might  well  be  influencing  the  flow  behavior  in 

these  tubes.   Since  no  drag  reduction  tests  were  made  using  the  larger 

tubes,  an  analysis  of  this  effect  was  not  performed. 

The  Midwest  differential  pressure  gage  used  to  measure  high 
test  section  differential  pressures  did  not  provide  the  precision 
or  range  required  in  this  study.   Its  continued  use  in  future  work  is 
not  recommended.   There  are  other  differential  pressure  gages  available 
which  will  provide  better  precision  for  a  given  range  of  pressure.   For 
future  work,  several  precise  differential  pressure  gages,  each  selected 
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U5 
to  cover  a  part  of  the  total  pressure  range  anticipated,  e.g.,  10,  50, 
and  100%  of  maximum  range,  should  be  used. 

It  is  recommended  that  a  means  of  controlling  the  temperature 
fluctuations  of  the  fluid  in  the  reservoir  be  provided.   These  fluctuations, 
caused  by  ambient  temperature  changes  in  the  laboratory,  made  it  necessary 
to  use  average  values  of  density  and  viscosity. 

2.  Performance  of  Vistanex  L-200  as  a_  Drag  Reducing  Additive   The 
results  of  this  study  indicate  that  a  marked  reduction  in  viscous  drag 
can  be  achieved  by  the  addition  of  minute  quantities  of  high  molecular 
weight  polyisobutylene  in  a  hydrocarbon  solvent.   Frictional  drag  de- 
creases of  as  much  as  65%  were  noted  when  turbulent  flow  experiments  were 
performed  with  a  50  wppm  concentration  of  Vistanex  L-200  dissolved  in  . 
military  jet  fuel.   The  friction  reduction  was  observed  to  increase  with 
increasing  flow  tube  diameter.   Further,  friction  reduction  increased 

with  increasing  Reynolds  numbers  up  to  a  certain  "critical"  Reynolds  number, 
after  which  friction  reduction  decreased  with  increasing  Reynolds  numbers. 
By  way  of  contrast ,  this  anomalous  behavior  was  not  observed  in  the 
laminar  region  of  fluid  flow  where  frictional  drag  was  slightly  higher 
than  that  predicted  by  the  Hagen-Poiseuille  relationship. 

3.  Sensitivity  of  Vistanex  L-200  to  Mechanical  Degradation  The 
solutions  of  polymer  tested  were  shown  to  be  extremely  sensitive  to 
mechanical  degradation.   Further,  repeated  testing  of  a  fluid  containing 

a  fixed  amount  of  polymer  indicated  that  the  drag  reducing  characteristics 
of  such  a  solution  could  be  irreversibly  destroyed.   The  degradation  was 
found  to  occur  when  the  solution  was  subjected  to  a  region  of  highly 
localized  shear  (a  needle  valve)  as  well  as  to  the  continuous  wall  shear 
in  the  experimental  tubes  employed. 

Any  practical  application  of  the  drag  reduction  phenomenon  must 
consider  the  limitations  imposed  by  degradation.   Therefore,  if  operational 
use  is  contemplated  in  a  full  scale  flow  system,  it  is  recommended  that 
mechanical  devices  producing  high  Reynolds  stresses  be  avoided  or,  alterna- 
tively, injection  of  the  additive  be  accomplished  downstream  from  such  a 
device. 
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4.  Relationship  Between  Diameter  Effect  and  Degradation  Although 
the  results  of  this  study  are  not  conclusive,  there  appears  to  be  a 
striking  similarity  between  what  has  been  termed  "the  diameter  effect" 
and  the  effects  of  mechanical  degradation.   It  is  therefore  recommended 
that  additional  work  be  undertaken  to  investigate  the  relationship 
between  the  flow  tube  diameter  effect  and  the  mechanical  degradation 
effect  of  drag  reduction. 

5.  Pseudo-Laminar  Flow  Characteristics  of  Fresh  Solutions  There 
were  strong  indications  in  the  experimental  work,  i.e.  the  pulsation  of 
the  pressure  measuring  devices,  that  a  suppression  of  turbulent  fluctu- 
ations might  well  be  one  of  the  major  characteristics  associated  with  the 
drag  reducing  influence  of  an  undegraded  polymer  solution.   However,  it 
is  recognized  that  the  pressure  measuring  equipment  employed  in  this 
study  could  not  permit  a  quantitative  verification  of  this  observation. 
It  is  therefore  recommended  that  additional  work  be  undertaken  to  verify 
these  qualitative  observations  with  pressure  sensors  with  high  frequency 
response,  e.g.,  strain  gages. 
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NOMENCLATURE 

Symbol  Definition 

d..  derivative  of  .. 

D  flow  tube  internal  diameter 

f  Fanning  friction  factor 

g  gravitational  constant 

K  correction  to  pressure  head  to  account  for  exit, 
entrance  and  kinetic  energy  losses 

L  flow  tube  test  section  length 

m  parameter  in  power  law  model 

n  index  of  non-Newtonian  behavior  in  power  law  model 

N_,  Deborah  number 
Db 

N  Reynolds  number 

P  pressure 

Q  volumetric  flow  rate 

R  radius  of  flow  tube 

V  average  velocity 

v  instantaneous  velocity 

x,y,z  cartesian  coordinates 

Greek  Symbols 

y  Newtonian  viscosity 

tt  3.1415... 

p  fluid  density 

x  shear  stress 

Tn  first  mode  relaxation  time  of  a  viscoelastic  solution 

Subscripts 

0,1,2...  evaluated  at  point  0,1,2... 

v  purely  viscous  fluid 

ve  viscoelastic  fluid 

w  evaluated  at  the  flow  tube  wall 

x,y,z  cartesian  coordinates 
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APPENDIX  A 

COMPUTER  PROGRAMS  USED  DURING  STUDY 

PROGRAM  NUMBER  1   COMPUTER  MODEL  OF  NEWTONIAN  FLUID  FLOW  IN 

SINGLE  PASS  VARIABLE  PRESSURE  FLOW  SYSTEM 

PROGRAM  NUMBER  2   COMPUTER  PROGRAM  TO  COMPUTE  CALIBRATION 

CONSTANT  OF  A  CAPILLARY  TUBE  VISCOMETER 

PROGRAM  NUMBER  3   COMPUTER  PROGRAM  TO  CALCULATE  INTERNAL 

'  DIAMETER  OF  FLOW  TUBES 


PROGRAM,  NUMBER  4   COMPUTER  PROGRAM  TO  ANALYZE  FLOW  DATA 
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COMPUTER  PROGRAM  NUMBER  1 

Computer  Model  of  Newtonian  Fluid  Flow  in  Single  Pass  Variable 
Pressure  Flow  System. 

This  program  was  developed  to  aid  in  the  operation  of  the  single 
pass  variable  pressure  flow  system  used  in  this  study.   It  can  also  be 
used  to  test  the  effect  that  variations  in  flow  tube  dimensions  such  as 
flow  tube  diameter,  test  section  length,  etc.,  will  have  on  the  performance 
of  the  flow  system.   A  sample  of  the  output  to  be  expected  from  this 
program  is  found  in  Table  2-1. 

Input  data  is  entered  in  seven  fields  each  containing  10  columns 
in  an  F10.0   format.   The  Reynolds  number  (RN)  for  which  performance 
characteristics  are  desired  is  entered  in  the  first  field.   If  the 
number  200  is  entered  in  this  field,  the  program  will  compute  and  write 
out  flow  characteristics  for  fluid  flow  at  various  Reynolds  numbers  in 
the  laminar  and  turbulent  flow  regimes  up  to  a  maximum  Reynolds  number 
of  140,000. 

The  flow  tube  diameter  (D)  in  inches  is  entered  in  the  second 
field.   Fluid  viscosity  (VISC)  in  centipoise  is  entered  in  the  third 
field.   Fluid  density  (SG)  in  grams  per  cubic  centimeter  is  entered  in 
the  fourth  field.   The  flow  tube  exit  length  (CLD),  entrance  length  (CLU) 
and  test  section  length  (L),  all  in  inches  are  entered  in  the  fifth,  sixth 
and  seventh  fields.   If  any  one  of  these  fields  is  left  blank,  the  program 
will  assume  an  entrance  length  of  75  diameters,  exit  length  of  25 
diameters,  or  test  section  length  of  100  diameters  for  the  blank  field. 

The  program  will  assist  in  operating  the  flow  system  at  a  given 
Reynolds  number  by  providing  the  minimum  weight  or  volume  of  fluid 
required  to  meet  the  criterion  that  the  flow  time  for  a  run  last  at 
least  100  seconds  and  by  providing  the  pressure  regulator  setting 
required  to  achieve  the  flow  regime  desired.   If  the  system  limitations 
of  available  pressure  or  reservoir  capacity  are  exceeded,  the  program 
will  automatically  signal  this  condition  to  the  operator.   A  blank 
card  at  the  end  of  the  input  data  will  exit  the  program  from  the  computer. 
Sample  input  data  are  located  at  the  end  of  the  program. 
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COMPUTER  PROGRAM  NUMBER  2_ 

Computer  Program  to  Compute  a  Calibration  Constant  for  a  Capillary 
Tube  Viscometer. 

This  program  was  developed  to  compute  an  average  calibration 
constant  for  a  Cannon-Fenske  Capillary  Tube  Viscometer  from  data  taken 
during  calibration  runs.   A  sample  of  output  to  be  expected  from  this 
program  is  found  in  Table  2-2. 

Input  data  is  entered  in  four  fields  each  containing  10  columns 
in  an  F10.0  format.   The  efflux  time  (TIME)  in  seconds  of  the  fluid 
flow  in  the  capillary  tube  is  entered  in  the  first  field.   The  fluid 
temperature  (TEMP)  in  degrees  Celcius  is  entered  in  the  second  field. 
The  fluid  density  (RHO)  in  grams  per  cubic  centimeter  and  fluid  vis- 
cosity (VIS)  in  centipoise  corresponding  to  the  measured  temperature 
are  entered  in  the  third  and  fourth  fields.   The  program  will  compute 
a  calibration  constant  for  each  calibration  run  and  the  arithmetic 
average  of  the  series  of  runs.   A  blank  card  at  the  end  of  the  input 
data  deck  will  exit  the  program  from  the  computer. 
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COMPUTER  PROGRAM  NUMBER  3_ 

Computer  Program  to  Calculate  the  Internal  Diameter  of  Flow  Tubes . 

This  program  was  developed  to  compute  an  average  internal  diameter 
for  the  flow  tubes  to  be  used  with  the  single  pass  variable  pressure  flow 
system  used  in  this  study.   A  sample  of  the  output  to  be  expected  from 
this  program  is  illustrated  in  Appendix  B. 

Each  series  of  data  cards  containing  information  taken  during 
flow  experiments  pertaining  to  a  tube  to  be  calibrated  is  preceded  by 
a  run  control  card.   The  run  control  card  contains  three  data  fields 
each  of  which  contains  10  columns  in  an  F10.0  format.   The  identifying 
number  (such  as ■ the  gage  number)  of  the  tube  is  entered  in  the  first 
field  (TUBE).   The  diameter  (SPEC)  of  the  tube  in  inches  as  supplied 
by  the  manufacturer  is  entered  in  the  second  field.   The  temperature 
(TEMP)  of  the  calibration  fluid  is  entered  in  the  third  field  in  degrees 
Celcius . 

The  data  cards  following  the  run  control  card  contain  information 
required  to  compute  the  internal  diameter  of  the  flow  tube  from  measure- 
ments taken  during  laminar  flow  runs.   Each  data  card  contains  seven 
data  fields  each  in  an  F10.0  format.   The  first  10  column  field  is 
left  blank.   The  test  section  length  (L)  in  inches  is  entered  in  the 
second  data  field.   The  elapsed  time  (T)  in  seconds  required  for  a 
predetermined  quantity  of  fluid  to  flow  through  the  tube  is  entered  in 
the  third  data  field.   The  value  of  fluid  density  (SG)  is  entered  in  the 
fourth  field  in  grams  per  cubic  centimeters.   The  value  of  fluid  vis- 
cosity (VISC)  is  entered  in  the  fifth  field.   This  value  may  be  either 
the  efflux  time  in  seconds  taken  from  the  calibrated  number  25/J680 
Cannon-Fenske  viscometer,  or  the  value  of  viscosity  in  centipoise  up  to 
a  maximum  of  5  cp.   The  volume  (VOL)  of  fluid  examined  during  the  flow 
experiment  expressed  in  milliliters  is  entered  in  the  sixth  data  field. 
The  differential  test  section  pressure  (P)  which  is  measured  in  inches 
of  CCL   is  entered  in  the  seventh  data  field.   Columns  71-80  are  left 
blank. 

A  blank  card  at  the  end  of  a  series  will  indicate  that  the  series 
has  been  completed  and  the  arithmetric  average  of  the  individually 


61 
computed  diameters  will  be  computed.   Two  blank  cards  at  the  end  of  a 

data  deck  will  exit  the  program  from  the  computer.   Sample  input  data 

are  located  at  the  conclusion  of  the  program  list. 
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COMPUTER  PROGRAM  NUMBER   U_ 

Computer  Program  to  Analyze  Flow  Data. 

This  program  was  developed  to  analyze  the  experimental  data  taken 
during  flow  experiments  with  the  variable  pressure  single  pass  flow 
system.   Each  series  of  data  cards  containing  information  taken  during 
flow  experiments  is  preceded  by  two  run  control  cards.   The  first  run 
control  card  contains  four  10  column  input  fields  in  an  F10.0  format. 

The  gage  number  of  the  flow  tube  used  during  the  flow  experiments 
is  entered  in  the  first  field  (TUBE).   The  density  in  grams  per  cubic 
centimeters  of  the  manometric  fluid  used  to  measure  differential  pressure 
is  entered  in  the  second  field  (SGMF).   The  flow  diameter  (D)  and 
test  section  length  (L)  both  in  inches  are  entered  in  the  third  and  fourth 
fields  respectively.   The  second  run  control  card  is  in  a  65A1  format 
on  which  a  general  description  of  the  experiment  is  entered. 

There  are  two  data  cards  required  to  record  information  pertinent 
to  each  flow  experiment.   The  first  data  card  contains  eight  10  column 
fields  in  an  F10.0  format.   The  first  two  fields  are  used  to  enter  the 
flow  tube  diameter  (D)  and  test  section  length  (L)  both  in  inches.   The 
third  field  is  used  to  enter  the  elapsed  time  (T)  of  the  flow  run  in 
seconds.   The  value  of  fluid  density  (SG)  is  entered  in  the  fourth  field 
in  either  CGS  units  up  to  5  grams  per  cubic  centimeter  or  British 
Engineering  System  units  for  values  exceeding  5  pounds  per  cubic  foot. 
The  value  of  fluid  viscosity  (VISC)  is  entered  in  the  fifth  field.   This 
value  may  be  either  the  efflux  time  in  seconds  taken  from  the  calibrated 
number  25/J680  Cannon-Fenske  viscometer,  the  temperature  in  °C  of 
distilled  water  if  used  in  the  flow  experiment  between  10-99°C  or  the 
value  of  viscosity  in  centipoise  of  the  fluid  if  less  than  10  cp.   The 
volume  (V0)  of  fluid  examined  during  the  flow  experiment  expressed  in 
milliliters  is  entered  in  the  sixth  data  field  if  a  volumetric  measure- 
ment is  taken.   The  differential  pressure  (P)  is  entered  in  the  seventh 
data  field.   This  value  may  be  either  in  inches  of  manometric  fluid  when 
the  V0  field  contains  a  digit,  or  pounds  per  square  inch  if  the  V0 
field  is  blank.   The  weight  (W)  in  pounds  of  fluid  examined  during  the 
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flow  experiment  if  a  weight  measurement  is  taken,  is  entered  in  the  eighth 

data  field. 

The  second  data  card  contains  three  10  column  fields  in  an  F10.0 
format.   It  is  used  to  enter  the  pressure  gage  calibration  correction 
(DEV)  and  maximum  (PLUS)  and  minimum  (SUB)  expected  error  bounds 
associated  with  the  pressure  measuring  instrument. 

A  blank  card  at  the  end  of  a  series  of  runs  will  indicate  that 
the  series  has  been  completed.  Two  blank  cards  will  exit  the  program 
from  the  computer. 

A  listing  of  data  taken  during  the  flow  experiments  is  located 
at  the  end  of  the  program  in  this  Appendix.   Examples  of  output  from 
this  program  are  found  in  Appendices  C  and  D. 
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APPENDIX  B 


COMPUTER  WRITE  OUT 
RESULTS  OE  THE  CALCULATIONS  TO  DETERMINE 
FLOW  TUBE  INTERNAL  DIAMETER 


8U 


CALIBRATION  OF  CAPILLARY  TUBE  GAGE  15. 
MANUFACTURERS  SUPPLIED  DIAMETER  IS  0.054 
PLUS  OR  MINUS  0;005  INCHES 


REYNOLDS  COMPUTED   MEASURED    FLOW     PRESSURE  OBSERVED   COMPUTED 

NUMBER    VISCOSITY  DENSITY     RATE     DIFFER.   FRICTION   DIAMETER 

(CP)     (GM/CC)   (GAL/MIN)  (IN  CCL4)FACT0R     (INCHES) 


80. 

20?. 

267. 

415. 

485. 

501  . 

680  . 

756. 

774. 
1007. 
10  81. 


0.9560 
0  .9560 
0.9560 
0  .9560 
0.9560 
0.9560 
0  .9560 
0.9560 
0 .9560 
0  .9560 
0.9560 


0.9975 
0 .9975 
0.9975 
0  .9975 
0  .9975 
0  .9975 
0.9975 
0.9975 
0.9975 

0  .9975 

0.9975 


0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


0.9 
2.3 

3.1 
4.7 
5.5 
5.6 
7.8 
8.6 
8.6 
11.3 
12.1 


ARITHMETIC  AVERAGE  0^  COMPUTED  INTERNA 


0.1992 
0  .0791 
0  .0599 
0.0385 
0.0330 
0.0320 
0.0235 
0.0212 
0.0207 
0.0159 
0.0148 
DI  AMETER 


0  .  0544 
0.0531 
0.0530 
0  .0532 
0 . 0530 
0.0536 
0  .053  0 
0.0531 
0.0534 
0.0533 
0.0533 
0.05331 
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CALIBRATION  OF  CAPILLARY  TUBE  GAGE  13. 
MANUFACTURERS  SUPPLIED  DIAMETER  IS  0.071 
PLUS  OR  MINUS  0.005  INCHES 


REYNOLDS  COMPUTED   MEASURED    FLOW     PRESSURE  OBSERVED   COMPUTED 

NUMBER    VISCOSITY  DENSITY     RATE     DIFFER.   FRICTION   DIAMETER 

(CP>     (GM/CC)   (GAL/MIN)  (IN  CCLOFACTOR     (INCHES) 


136. 

255. 

362. 

512. 

700. 

870. 
1002. 
1257. 
1396. 
1598. 
1800. 
1915. 
2055. 


0 .9840 
0.9840 
0 .9840 
0 .9840 
0.9840 
0 .9840 
0 .9840 
0.9840 
0.9840 
0 .9840 
0.9840 
0 .9840 
0.9840 


0.9973 

0 ,9973 

0 .9973 

0 .9973 
0 .9973 
0.9973 
0 .9973 
0 .9973 
0 .9973 
0 .9973 

0 .9973 
0 .9973 
0  .9973 


0.0 

0.0 

0.0 

0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


0.8 

1.5 

2.3 

3.5 

4.5 

5.4 

6.5 

8.2 

9.0 

10.2 

12.1 

13.0 

14.1 


0.1176 
0.0628 
0  .  0442 
0.0312 
0 .0228 
0 .0184 
0 .0160 
0.0127 
0.0115 
0.0100 
0.0089 
0.008^ 
0.0076 


0.0733 
0  .0732 
0  .0717 
0.0702 
0.0718 
0,0725 
0.0716 
0  .0713 
0.0715 
0.0720 
0.07Q7 

0.07Q5 
0.07Q2 


ARITHMETIC  AVERAGE  OE  COMPUTED  INTERNAL  DIAMETER   0.07158 
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CALIBRATION  OF  CAPILLARY  TUBE  GAGE  12. 
MANUFACTURERS  SUPPLIED  DIAMETER  IS  0.085 
PLUS  OR  MINUS  0.005  INCHES 


REYNOLDS  COMPUTED   MEASURED    FLOW     PRESSURE  OBSERVED   COMPUTED 

NUMBER    VISCOSITY  DENSITY     RATE     DIFFER.   FRICTION   DIAMETER 

CCP)     (GM/CC)   (GAL/MIN)  (IN  CCL4)FACT0R     (INCHES) 


178. 

0 .9828 

0 

.9972 

0.0 

0.9 

0 .0898 

0.0834 

604. 

0 .9828 

0 

.9972 

0.0 

2.5 

0.0265 

0.0878 

971. 

0 .9826 

0 

.9972 

0.0 

4  .  3 

0.0165 

0 .0865 

1225. 

0.9828 

0 

.9972 

0.0 

5.5 

0.0131 

0  .0859 

1427. 

0 .9828 

0 

.9972 

0.0 

6.3 

0.0112 

0.0863 

1621. 

0.9826 

0  , 

.9972 

0.0 

7.3 

0.0099 

0,0856 

1815. 

0 .9828 

0 

,9972 

0.0 

8.2 

0.0088 

0.0854 

2014. 

0 .9628 

0 

.9972 

0.1 

9.2 

0.0079 

0  .0852 

2149. 

0 .9628 

0 

.9972 

0.1 

10.1 

0.0074 

0  .0844 

2154. 

0 .9626 

0 

.9972 

0  .1 

10.1 

0.0074 

0  .0844 

ARITHMETIC  AVERAGE  Of    COMPUTED  INTERNAL  DIAMETER   0.08548 
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CALIBRATION  OF  CAPILLARY  TUB-P  GAGE  11, 
MANUFACTURERS  SUPPLIED  DIAMET-R  IS  0.094 
PLUS  OR  MINUS  0.005  INCHES 


REYNOLDS  COMPUTED   MEASURED    PLOW     PRESSURE  OBSERVED   COMPUTED 

NUMBER    VISCOSITY  DENSITY     RATE     DIFFFR.   FRICTION   DIAMETER 

(CP)     (GM/CC)   (GAL/MIN)  (IN  CCL4)FACT0R     (INCHES) 


455. 

639. 

950. 
1173. 
1412. 
1667. 
1896. 
2029. 


0.9352 
0.9352 
0.9352 
0.9352 
0.9352 
0.9352 
0.9352 
0 .9352 


0.9974 

0 .9974 

0 .9974 
0 .9974 
0  .9974 
0 .9974 
0 . 99/4 

0 ,9974 


0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.1 
0.1 


1.5 
2.2 
3.4 
4.3 
5.1 
6.1 
7.0 
7.7 


0.0352 
0.0250 
0.0168 
0.0136 
0  .  0  1 1  3 


0.0953 
0 , 0939 
0.0927 
0  .0920 
0  ,0924 


0.0096  0.0920 
0.0084  0.0917 
0.0079    0.0911 


ARITHMETIC  AVERAGE  Of  COMPUTED  INTFRNAL  DIAMETER   0.09?64 


CALIBRATION1  OF  CAPILLARY  TUBE  GAGE  10. 
MANUFACTURERS  SUPPLIED  DIAMETER  IS  0.106 
PLUS  OR  MINUS  0.005  INCHES 


REYNOLDS  COMPUTED   MEASURED    FLOW     PRESSURE  OBSERVED   COMPUTFD 

NUMBER    VISCOSITY  DENSITY     RATE     DIFFER.   FRICTION   DIAMETER 

fCP)  (GM/CO       (GAL/MIN)     (IN    CCL4)FACT0R  (INCHES) 


358. 

0.9351 

0 

.9973 

0.0 

0.9 

0.0447 

0,1084 

656. 

0.9363 

0 

.9985 

0.0 

1.8 

0.0244 

0.1063 

817. 

0.9351 

0 

.9973 

0.0 

2.1 

0.0196 

0.1076 

947. 

0.9351 

0 

,9973 

0.0 

2.5 

0.0369 

0.1066 

1164. 

0.9363 

0 

.9985 

0.0 

2.9 

0.0138 

0.1082 

1309. 

0.9351 

0 

,9973 

0.0 

3.6 

0.0122 

0 .1051 

1558. 

0.9351 

0 

.9973 

0.0 

4.4 

0.0103 

0.1042 

1833. 

0.9351 

0 

.9973 

0.1 

5.3 

0.0087 

0.1034 

2034. 

0.9351 

0 

.9973 

0.1 

5.9 

0.0079 

0.1033 

ARITHMETIC  AVERAGE  OF"  COMPUTED  INT-RNAL  DIAMETER   0.10590 
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APPENDIX  C 


COMPUTER  WRITE  OUT 
RESULTS  OF  FLOW  DATA  TAKEN 
TO  DETERMINE  INTERNAL  CHARACTERISTICS  OF  FLOW  TUBES 


90 
ANALYSIS  OF  FLOW  DATA 

flow  of  tap  water  to  determine  intfrior  roughness  characteristics 

in  flow  tube  number  i 5 . 


TUBE  DIMENSIONS  (INCHPS)   FLOW  TUBE  DIAMETER  0.0533 

TEST  SECTION  LENGT^   5.406 


REYNOLDS  COMPUTED   MEASURED    FLOW1 
NUMBER    VISCOSITY  DENSITY     RATE 

(Co.)'     (6M/CC)   (GAL/MIN) 


PRESSURE  OBSERVED   COMMUTED 
DIFFFP,   FRICTION   FRICTION 
(PSD    FACTOR     FACTOR 


65350.     0.976? 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 

64731.     1.0718 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0.997?       1 

DUE  TO  U N C E R T  A 
DLE  TO  UNCERTA 

0  .  9  9  h  7  1 

DUE  TO  UN CERT  A 
DUE  TO  UNCERTA 


53122.     1.0  773    0.996  7       0 

MAX  EXPECTED  VALUES  DLE  TO  UNCERTA 
MIN  EXPECTED  VALUES  DLE  TO  UNCERTA 


42013.     1.0773 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0.9987  0 
DUE  TO  UNCERTA 
DLE  TO  UNCERTA 


29341  .  1.0823  0  .9987  0 
MAX  EXPECTED  VALUES  DUE  TO  UNCERTA 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTA 

23162.  1.0828  0.9987  0 
MAX  EXPECTED  VALUES  DUE  TO  UNCERTA 
MIN  EXPECTED  VALUES  DUE  TO  U N C E R'T A 

19  5  0  4.  1.0800  0.9967  0 
MAX  EXPECTED  VALUES  DLE  TO  UNC=PTA 
MIN  EXPECTED  VALUES  DUE  TO  U N C E R T A 


08 

NTY 

NTY 

17 

NTY 
NTY 

97 

NTY 
NT  V 

76 
NT  v 
NTY 

54 

NTY 

NTY 

4? 

NTY 

NTY 

3  6 

\t  J  Y 
n;TV 


362 
364 

356 

413 
41  6 

40  6 


100 
101 

98 

66 
67 
65 


298.4 
290.2 
285  .  6 

169.6 
191.3 
167  .2 


51.5 

51  .9 
50.5 


0 . 00554 
0 .00558 

0.00545 

0  .  00536 
0  ,00539 
0  .00529 

0.00549 
0.0055? 
0  .00543 

0  .00577 
0.0058? 
0,00570 

0  .  0062O 
0  .  00626 
0  ,0060* 

0.00663 
0  .  00669 
0.00651 

0.00723 
0.00729 

0  .00709 


0  .00495 


0.00496 


0.00521 


0  .00552 


0.00604 


0 . 00641 


0 .00669 
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analysis  of"  plow  data 

flow  of  pure  jp4  to  determine  interior  roughness  characteristics 

in  flow  tube  number  15. 


FLOW  TUBE  DIMENSIONS  (INCHES)   FLOW  TURP  DIAMETER  0.0533 

TEST  SECTION  LENGTH   5.406 


REYNOLDS  COMPUTED   MEASURED    FLOW     PRESSURE  OBSERVED   COMPUTED 
NUMBER    VISCOSITY  DENSITY     RATE     DIFFER,   FRICTION   FRICTION 
(CD)     CGM/CC)   (6AL/MIN)   (PSI)    FACTOR     FACTOR 


4  4  2  5  7  .     0.7433    0.7  64  4       0 

MAX  EXPECTED  VALUES  DUE  TO  UwCERJA 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTA 


37735  . 
MAX  EXPE 
HI N  EXPECTED 


0.7  4  33    0.7  64  4       0 

C T E D  VALUES  DUE  TO  U N C E R T  A 

VALUES  DUE  TO  UNCERTA 


?9464.  0.7433  0.7644  0 
MAX  EXPECTED  VALUES  DUE  TO  U N C E R T A 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTA 


29270.     0.794? 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 

254  8  9.     0.7942 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 

22118.     0.794? 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


16120  . 
MAX  EXPE DTE 
MIN  EXPPCTE 


0.7942 
[3  V  A  L  U  u  S 
D  VALUES 


0.7657  0 
DUE  TO  UNCERTA 
DUE  TO  UNCERTA 

0.7657       0 

DUE  TO  UNCERTA 
DUE  TO  UNCERTA 

0.  7657       0 

DUE  TO  UNCERTA 
DUE  TO  UNCERTA 

0 .7657  0 
CUE  TO  UNCERTA 
DUE  TO  UNCERTA 


73 

N  T  Y 

NTY 

62 
NTY 

NTY 

48 

MTY 

NTY 

51 
NJY 

NTY 

45 
MTY 

N  T  Y 

39 

NTY 
\TY 

28 

NTY 
NTY 


128 

12° 
125 

97 
98 
95 


69 
70 
67 


63.  0 
63.6 
61.6 


55.7 
56.  4 
54 .0 

41  .4 
41.8 
40.2 

21  .  3 
22.0 
21  .2 


0  .00565 
0.00569 
0 . 00554 

0.00592 

0.00597 
0.00531 

0.00626 
C ,0063? 
0.00613 

0.00611 
0.00619 
0 . 0059? 

0 .00649 
0 . 00658 
0.00630 

0.00641 
0 . 00643 
0.0082? 

0 .00635 
0  .  0064? 
0.00617 


0.00545 


0.00568 


0.00604 


0.00605 


0.006?6 


0.00649 


0.00702 
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ANALYSIS  OF  FLOW  DATA 

FLOW  OF  PURE  Jp4  TO  DETERMINE  INTERIOR  ROUGHNESS  CHARACTERISTICS 

IN  F|_Oa'  TJBc  NUMBER  15. 


FLOW  TURE  DIMENSIONS  (INCHES) 


FLOW  TUBE  DIAMETER  0.0533 
TEST  SECTION  LENGTH   5,406 


REYNOLDS  COMPUTED   MEASURED    FLOW     PRESSURE  OBSERVED   COMPUTED 
NUMBER    VISCOSITY  DENSITY     RATE     DIFFER.   FRICTION   FRICTION 
( C  p > .     ( G  M / C  C 1   ( G  A  L / M I N )   C  PS  M    FACTOR     FACTOR 


39707.     0.7295 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0.7613  0 
DUE  TC  UNCERTA 
DUE  TO  UNCERTA 


5371?.  0.7295  0.7613  0 
MAX  EXPECTED  VALUES  DUE  TO  UNCERTA 
MIN  EXPECTED  VALUES  DUE  TO  UNCERJA 


73536.     0.7295 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 

86633.     0.7295 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 

62431.     0.7*33 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0.7613       1 

DUE  TO  UNCERTA 
DUE  TO  UNCERTA 

0.7613       1 

DUE  TO  UNCEPTA 
DUE  TO  UNCERTA 

0.7644  1 
DUE  TO  UNCEPTA 
DUE  TO  UNCERTA 


5337?.     0.7  4  33    0.7644       0 

MAX  EXPECTED  VALUES  DUE  TO  UNCERTA 
MIN  EXPECTED  VALUES  DUE  TC  UNCERTA 


64 
\!  T  Y 
NTY 

87 

NTY 

NTY 

19 

NTY 

NTY 

40 

NiTY 

NTY 

02 

VTY 
NTY 

88 
MTY 

VTY 


102.5 
10  3,3 
10  0.1 

176.8 
173.9 


314 
317 
310 

437 
43^ 
47? 


247.4 
250.4 
244  .  0 

134  .3 

136.1 
131.2 


0.0  0  530   0.00560 

0.00583 

0.00567 

0.0  0  547   0.00520 

0.0055* 

0.00539 


0  0  519   0.00480 
00524 


0. 

0. 
0.00512 


00520   0.00461 
0  0  5  2  3 


0. 
0. 

0.00514 


0.00547   0.01500 
0.00554 

0.00540 


0  0  5  5  8   0.00520 
00564 


0. 

0  . 

0  .  0054Q 
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analysis  of  plow  data 

flow  op  tap  water  to  determine  interior  roughness  characteristics 

in  flow  t u h e  n u m per  i 3 . 


FLOW  TURE  DIMENSIONS  (INCHES) 


FLOW  TUBE  DIAMETER  0.0716 
TEST  SECTION  LENGTH   7,109 


REYNOLDS  COMPUTED   MEASURED 
NUMBER    VISCOSITY  DENSITY 
(  C  P. )     (  G  M  /  C  C  ) 


FLOW     PRESSURE  OBSERVED   COMPUTED 
RATE     DIFFF9.   FRICTION   FRICTION 
(GAL/MIM)   <  PS  I  )    FACTOR     FACTOR 


966 Sn.     1.0664 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0.9967  2.3  4  434.5  0  .  0  0  4  6  R 
DUE  TO  UNCERTAINTY  437.  4  0.0.0472 
DUE  TO  UNCERTAINTY    430.9    0.00465 


0  .00449 


83247.  1.0691  0,9967  2.02  330.7  0.00479 
MAX  EXPECTED  VALUES  DUE  TO  UNCERTAINTY  332.6  0,0  0482 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTAINTY    327.0    0.00474 


0.00466 


72868.     1.06  91 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0.99  6  7  1.77  258.0  0.00488 
DUE  TO  UNCERTAINTY  259.6  0.00491 
DUE  TO  UNCERTAINTY    255.7    0.0 048 3 


0  .004*1 


61962.  0.9656  0.996Q  1.39  162.7  0.00500 
MAX  EXPECTED  VALUES  DUE  TO  UNCERTAINTY  164.2  0.00505 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTAINTY    160.6    0. 00494 


0  .00501 
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ANALYSIS  OF  FLOW  DATA 

flow  of  pure  jj4  to  determine  interior  roughness  characteristics 

in  flow  tube  number  i?. 


FLCH  TUBE  DIMENSIONS  (INCHES) 


FLOW  TUB-  DIAMETER  0.0716 
TEST  SECTION  LENGTH   7.10? 


REYNOLDS  COMPUTED   MEASURED    FLOW 
NUMBER    VISCOSITY  DENSITY     RATE 

(CP)     (GM/CC)   (GAL/MIN) 


PRESSURE  OBSERVED   COMPUTED 
DIFFER.   FRICTION   FRICTION 
( P  S  I  )    fACTOR     FACTOR 


9980  0.     0.7124    0.7610       2.12 
MAX  EXPECTED  V A l U & S  D u E  rQ  UNCERTAINTY 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTAINTY 


74795.  0.7124  0.7610  1 
MAX  EXPECTED  VALUES  DUE  TO  UNCERTA 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTA 


50396.     0.7124 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0.7610  1 
DUE  TO  UN CERT  A 
DUE  TO  UNCERTA 


55823.     0.7646  0.7644       1 

MAX  EXPECTED  VALUES  DUE  TO  UN CERT A 

MIN  EXPECTED  VALUES  DUE  TO  UNCERTA 

42670.     0.7646  0.7644       0 

MAX  EXPECTED  VALUES  DUE  TO  UNCERTA 

MIN  EXPECTED  VALUES  DUE  TO  UNCERTA 


31227.     0.764  6 

MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


37016.     0.764 
MAX  EXPECTED  VA 
MIN  EXPECTED  VALUES 


UES 


0 , 7644  0 
DUE  TO  UNCERTA 

DUE  TO  UNCERTA 

0 . 7644  0 
DUE  TO  UNCERTA 
DUE  TO  J N C E R T A 


26231.  0.7646  0.7644  0 
SAX  EXPECTED  VALUES  DUE  TO  UNCERTA 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTA 


59 
N  T  Y 
NTY 

07 
N  T  Y 
NTY 

26 

\>TY 

NTY 

97 

VTY 
NTY 

71 

N  T  Y 

NTY 

84 

NTY 

NTY 

59 

vTY 
v  T  Y 


269.3 
272.6 
26*.  7 

15  9.1 
160.7 
156.0 


77 

78 
75 

107 
10  3 
105 


67.0 
67,9 
65.0 


37.3 
37.7 

36.1 

51.5 

50.0 

25.9 
26.2 
25.2 


0  .00466 
0,00473 

0 .00461 

0.00490 
0.00495 
0 .00480 

0 . 00523 

0.00529 
0.00509 

0  .00519 
0 . 00523 
0.00506 

0 . 0055^ 
0.00560 
0.00536 

0.00574 
0 .00580 
0.00557 

0.0  0565 

0.0057? 
0.00548 

0 .0056- 
0.0057? 
0 .0055H 


0.00445 


0.00478 


0 . 00528 


0 .00515 


0.00550 


0.00595 


0.00570 


0.00622 
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ANALYSIS  OP  FLOW  DATA 
FLOW  OF  TAP  WATER  TO  DETERMINE  INTERIOR  ROUGHNESS  CHARACTERISTICS 

in  flow  tube  number  m. 


FLOW  TURE  DIMENSIONS  (INCHES) 


FLOW  TUBE  DIAMETER  0  .  10  5  9 

TEST  SECTION  LENGTH  10.609 


REYNOLDS  COMPUTED   MEASURED 
NUMBER    VISCnSITY  DENSITY 
(CD->     (G^'/CC) 


FLOW 
RATE 
(  G  A  L  /  M  I  N  ) 


PRESSURE  OBSERVED   COMPUTED 
DIFFER.   FRICTION   FRICTION 
( P S  I  )    FACTOR     FACTOR 


7571?.     1.022?    0.99  7  5       2.60 
MAX  EXPECTED  VALUES  DUE  TO  UNCERTAINTY 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTAINTY 


124 
125 
122 


0  .  00518 
0 . 00522 
0  .00507 


0.00477 


5  6  710.     1.0147 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0.99  7  5  1.93  76.0  0.00570 
DUE  TO  UNCERTAINTY  76.7  0.00576 
DUE  TO  UNCERTAINTY     74,1    0.00*57 


0.00513 


4547  7.     1.0172    0.9975       1.55 
MAX  EXPECTED  VALUES  DUE  TO  UNCERTAINTY 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTAINTY 


5?  .  5 
52  .  9 

51.1 


0  .  00610 
0  .00615 
0.00593 


0.0054? 


115435.     1.0638 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0.9975       4.12 
DUE  TO  UNCERTAINTY 
DUE  TO  UNCERTAINTY 


312 
31^ 

309 


0.0051* 
0  .00518 
C  .  0  0  5 1 0 


0  .  004?9 


84734.  1.0534  0.9975  3.01  170.5  0.00527 
MAX  EXPECTED  VALUES  DOE  TO  UNCERTAINTY  171.9  0.00531 
MIN  EXPECTED  VALUES  DUE  TC  UNCERTAINTY    168.2    0.00520 


0  .00464 
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ANALYSIS    OF    FlOW    DATA 
FLOW    Of    PURE    JP<    TO    DETERMINE    I.NTERIOR    ROUGHNESS    CHARACTERISTICS 

in  flow  tube  number  i o , 


FLOW  TUBE  DIMENSIONS  (INCHES) 


FLOW  TURF  DIAMETER  0.1059 
TEST  SECTION  LENGTH  10.609 


REYNOLDS  COMPUTED   MEASURED    FLOW 
NUMBER    VISCOSITY  DENSITY     RATE 

(CD)     (GM/CC)   (GAL/MIM) 


PRESSURE  OBSERVED   COMPUTED 
DIFFcp.   FRICTION   FRICTION 
(PS  I)    FACTOR     FACTOR 


90630.  0.7*75  0.7633  2 
MAX  EXPECTED  VALUES  DuE  TO  UNCERTA 
MIN  EXPECTED  VALUES  DUE  TC  UNCERTA 


73626.     0.7475 

MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0.7633       2 

DUE  TO  U\CERTA 
DDE  TO  UNCERTA 


56591.     0.7647  0.7638       1 

MAX  EXPECTED  VALUES  DUE  TO  UNCERTA 

MIN  EXPECTED  VALUES  DUE  TO  UNCERTA 

46015.     0.7647  0.7635       1 

MAX  EXPECTED  VALUES  CUE  TO  UNCERTA 

MIN  EXPECTED  VALUES  DUE  TO  UNCERTA 


340&5.     0.7647 
MAX  EXPECTED  VALUE'S 
MIN  EXPECTED  VALUES 


0.7636       1 

DUE  TO  UNjCEPTA 
DUE  TO  UNCERTA 


98 
NTY 

NTY 


*2 

N:TY 

NTY 

90 

NTY 

NTY 

54 

N:  T  V 

NTY 

14 
N  T  Y 

NTY 


128 

126 

87 
68 
65 

6  3 

64 
62 

3  8 
38 
37 


20  .3 
21.0 
20.2 


0  .00532 
0,00538 

0.00520 

0.00548 
0  .  00555 
0,00534 

0  .00650 
0 . 0065° 
0.00630 

0 .00589 
0  .  0  Q  5  9  5 
0  .00571 

0  .00582 
0  .00588 
0 .00565 


0.00456 


0.0O4R0 


0.00513 


0.005*0 


0  .00562 
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ANALYSIS  OF  FLOW  DATA 
A  M  I  N  A  R  F  L  0  a'  TEST  USING  P  U  »  E  T  A  p  WATER 

in  flow  tube  number  13, 


FLOW  TUBE  DIM-ENS  IONS  (INCHES)   FLOW  TUB?  DIAMETER  0.0716 

TEST  SECTION  LENGTH   7.109 


REYNOLDS  COMPUTED   MEASURED    FLOW     PRESSURE  OBSERVED   COMPUTED 
NUMBER    VISCOSITY  DENSITY     RATE     DIFFEP.   FRICTION   FRICTION 
(  C  P  )     (  G  M  /  C  C  )   (  G  A  L  /  M  I  N  )   (PSD    FACTOR     FACTOR 


36?  . 

502  . 

703. 

881. 
1002. 
1252. 
1395. 
1607. 
17  7  6. 
1864  . 
2016. 


0  .  9840 
0  .9840 
0.9840 
0  .9840 
0 .9640 
0  .9840 
0  .9840 
0  .  984C 
0  .98^0 
0  .  9840 
0.9640 


0  .9973 
0.9973 
0.9973 

0  .9973 

0  .  99/3 
0.9973 

0 .9973 
0.9973 
0.9973 
0 .9973 
0 .9973 


0.01 

0.01 
0  .02 
0,02 
0.02 
0  .03 
0  .03 
0.04 
0  .04 
0.0  4 
0.05 


0  .0 
0.1 
0.1 
0  .1 
0.1 
0.2 
0.2 
0  .  2 
0.3 
0.3 
0.3 


0.04398 
0.03443 
0  .  02245 
0  .01730 
0  .01598 
0.01298 
0.01151 
0  .0097* 
0  .00947 
0  .00°05 
0 . 00657 


0.04418 
0.03185 
0. 02276 
0  .  0181.6 
0  .  01597 
0.01278 
0  .01147 
0.00996 
0  .00900 
0.00849 
0.00794 
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ANALYSIS  OF  FLOW  DATA 
LAMINAR  FLOW  TEST  USING  PURE  Tap  WATER 

in  flow  tube  n u m b e r  is, 


FLO1-'  TUBE  DIMENSIONS  (INCHES)   FLOW  TURF  DIAMETER  0.0533 

TEST  SECTION  LENGTH   5.406 


REYNOLDS  COMPUTED   MEASURED    FLOW     PRESSURE  OBSERVED   COMPUTED 
NUMBER    VISCOSITY  DENSITY     RATE     DIFFER.   FRICTION   PRICTION 
(CP)     (GM/CC)   (GAL/MIN)   (PSI)    FACTOR     FACTOR 


265. 

0 

,9560 

0.9975 

0.00 

0.1 

414  , 

0 

.9560 

0.9975 

0.01 

0.1 

676. 

0 

.9560 

0 .9975 

0.01 

0  .? 

482. 

0 

.9560 

0.9975 

0.01 

0.1 

503. 

0 

.  9  5  6  0 

0.99  7  5 

0 .01 

0.1 

775. 

0 

.9560 

0.9975 

0 .01 

0.2 

75?. 

0 

.9560 

0 .9975 

0,01 

0.2 

1081  . 

0 

.9560 

0 .9975 

0  .02 

0.3 

1007. 

0 

.9560 

0.9975 

0.02 

0  .2 

0.06172  0.0  6029 

0.03907  0.03864 

0,0242*  0.02368 

0.0 3 390  0.03318 

0.0 31 2*  0.03182 

0.02050  0.02063 

0.02166  0.02127 

0.01433  0.01461 

0.01587  0.01583 
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APPENDIX  D 


COMPUTER  WRITE  OUT 
RESULTS  OE  FLOW  DATA  TAKEN 
DURING  DRAG  REDUCTION  EXPERIMENTS 
(FLOW  OF  50  WPPM  PIB  IN  MILITARY  UET  FUEL) 
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analysis  of  flow  data 

flow  of  jet  fucl  with  degraded  pit  (fourth  pass) 
in  flow  tube  mumper  15. 


FLOW  TUPE  DIMENSIONS  (INCHES) 


FLOW  TUB^:  DIAMETER  0.0533 
TEST  SECTION  LENGTH   5.406 


REYNOLDS  COMPUTED   MEASURED    FLOW 
NUMBER    VISCOSITY  DENSITY     PATE 

(CP)-    (GM/CC)   (GAL/MIN) 


PRESSURE  OBSERVED   COMPUTED 
DIFFER.   FRICTION   FRICTION 
(PSp    FACTOR     FACTOR 


80079.     0.77  3  9 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0,7685       1.36     319,6    0.00399 
DUE  TO  UNCERTAINTY    32  3.4    0.00*04 

DUE  TO  UNCERTAINTY    316.4    0.00395 


0.00470 


562  76.     0.7  73  9 

MAX  EXPECTED  VALUES 
MIN  EX=ECTEn  VALUES 


0.7685  0.96  156.1  0.00395 
DUE  TO  UNCERTAINTY  157,6  0.00398 
DUE  TO  UNCERTAINTY    153.0    0.00387 


0.00514 


43900.     0.7739 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0.7685       0,75 
DUE  TO  UNCERTAINTY 
DUE  TO  UNCERTAINTY 


9  °  .  2 
1 0  n  .  2 

97  .0 


0.00412 
0.0041 6 

0.00403 


0 . 0P546 
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ANALYSIS  OF  FLOW  DATA 
FLOW  0^"  JET  FUEL  WITH  DEGRADED  PI8  (THIRD  PASS) 

in  flow  tube  number  i5, 


FLOW    TIJRE    DIMENSIONS     (INCHES) 


FLOW  TUBE  DIAMETER  0  .  0  5)33 
TEST  SECTION  LENGTH   5.406 


REYNOLDS  COMPUTED   MEASURED 
NUMBER    VISCOSITY  DENSITY 
(CP)     (GM/CC) 


F  L  0  W 
RATE 
(  G  A  L  /  M  I  N  ) 


PRESSURE  OBSERVED   COMPUTED 
DIEFFP.   FRICTION   FRICTION 
(PSD    FACTOR     FACTOR 


78847.     0.7654 

MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0.7672       1.33 
DUE  TO  UNCERTAINTY 

DUE  TO  UNCERTAINTY 


2Q-1 
278 


0.00370 
0.00374 

0  .003  66 


0.00472 


60353.     0.7*54 
MAX  EXPCCTED  VALUES 
MIN  EXPECTED  VALUES 


0,767p  1.02  163.2  0.00  3  66 
DUE  TO  UNCERTAINTY  164.8  0.00370 
DUE  TO  UNCERTAINTY    16  0.0    0.00359 


0.00505 


47755.  0.7  6  54  0.767?  0.80  101.3  0.00363 
MAX  EXPECTED  VALUES  DUE  TO  UNCERTAINTY  102.3  0.003*7 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTAINTY     99.0    0,00355 


0.00535 


102 


analysis  of  n.ow  data 

flow  or  jet  fuel  with  degraded  fir  (second  pass) 
in  flow  tube  number  ip, 


FLOW  TU^E  DIMENSIONS  (INCHPS)   FLOW  TURF  DIAMETER  0.0533 

TEST  SECTIOM  LENGTH   5.406 


REYNOLDS  COMPUTED   MEASURED    FLOW     PRESSURE  OBSERVED   COMPUTED 
NUMBER    VISCOSITY  DENSITY     RATE     DIFFER.   FRICTION   FRICTION 
( C P )  ■    ( G M / C C  J   ( G A L / M I N )   ( P S I )    FACTOR     FACTOR 


8479. 

0.7551 

0 .7  569 

0.14 

5.6 

12084  . 

0.7531 

0,7669 

0  .20 

9.1 

14941  . 

0  .7581 

0 .7669 

0  .25 

12.7 

17926. 

0.7531 

0 , 7669 

0.30 

16.1 

2  0  212. 

0 .7581 

0 .7669 

0.34 

IP. 5 

22402  . 

0  .7531 

0  .7669 

0  .37 

2  2.7 

0.00647  0.O0B24 

0.00521  0.00754 

0.00473  0.00715 

0.00419  0.00694 

0.00377  0.0  0663 

0.00377  0.00647 
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ANALYSIS  OF  Ft  OW  PAT  A 

FLOW  OF'  5  0  W  P D  M  Pig  IN  JET  FUEL  (FIRST  PASS  THRU  ORIGINAL  SYSTEM) 

in  flow  tube  number  is. 


FLOW  TUBE  DIM: MS  I  CMS  (INCHES)   FLOW  TUBE  DIAMETER  0.0533 

TEST  SECTION!  LENGTH   5,405 


REYNOLDS  COMPUTED   MEASURED    FLOW     PPESSURR  OBSERVED   COMPUTED 
NUMBER    VISCOSITY  DENSITY     RATE     DIFFER,   FRICTION   FRICTION 
(CP).    (GM/CC1   (GAL/MIN)   <  P  ?  I  )    FACTOR     FACTOR 


25802. 

0.761? 

0  .7  65? 

0.^3 

19.4 

18314  . 

0.761? 

0.7652 

0  .31 

13  .5 

11351  . 

0.7612 

0  .765? 

0.19 

5  .6 

12875. 

0.761? 

0 .765? 

0.22 

6.8 

28379. 

0.761? 

0 . 76?2 

C.^8 

23  ,9 

0.00240  0.00  624 

0.00?  8 3  0.00680 

0.00360  0.007  6  6 

0.00336  0.00743 

0  .  0  0  2  ?  '>  0.0  0609 
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ANALYSIS  OF  FLOW  n a T A 

FLOW  0 F  5 0 W P P M  PIS  IN  J E T  FUEl(FIRST  pas?  THRU  ORIGINAL  SYSTEM) 

I N  F L 0 W  TUBE  NUMBER  15, 


FLOW  TUPE  DIMENSIONS  (INCHES) 


FLOW  TUBE  DIAMETER  0.0533 
TEST  SECTION  LENGTH   5.406 


REYNOLDS  COMPUTED   MEASURED    FLOW 
NUM8ER    VISCOSITY  DENSITY     RATE 

(CP).    (GM/CC)   (GAL/MIN) 


PRESSURE  OBSERVED   COMPUTED 
DIFFED.   FRICTION   FRICTION 
( P S  P    FACTOR     FACTOR 


5300  5.  0.7  7  05  0.7670  0.90  78.4  0.00225 
MAX  EXACTED  VALUES  DUE  TO  UNCERTAINTY  79.3  0.00227 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTAINTY     76.0    0.00213 


0  .  0T521 


42242.     0.7705    0.76  7  0       0.72 
MAX  EXPECTEO  VALUES  DUE  TO  UNCERTAINTY 

MIN  EXPECTED  VALUES  DUE  TO  UNCERTAINTY 


40  .3 
40  ,9 

39.0 


0.00182 
0.00185 
0  .00176 


0 .00552 


33057.  0.7705  0.767Q  0,56  2K.9  0.00191 
MAX  EXPECTED  VALUES  DUE  TO  UNCERTAINTY  26,2  0.00194 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTAINTY     25.0    0.00184 


0.00587 


30173.     0.7705 

MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0  .  767Q       0  ,51      20.7    0  .00184 
DUE  TO  UNCERTAINTY     21.0    0.00186 

DUE  TO  UNCERTAINTY     2^.0    0.00177 


0.00600 


23433.     0.7705 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0.76  7  0  0.^0  13.0  0.00139 
DUE  TO  UNCERTAINTY  13.1  0.00192 
DUE  TO  UNCERTAINTY     12.5    0.00183 


0.00639 
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ANALYS  is  of  fi  ow  data 

FLOW  OF  5 0 W P P M  P I B  IN  J E T  F U F L ( F  I  R S T  PASS  THRU  ORIGINAL  S Y S T e M  ) 

in  flow  tube  number  i  3 , 


FLO'-1  TURE  DIMENSIONS  (INCHES) 


FLOW  TUBE  DIAMETER  D.Q716 
TEST  SECTION  LENGTH   7.109 


REYNOLDS 

NUMBER 


COMPUTED 
VISCOSITY 
(CP)  . 


MEASURED 

DENSITY 

(GM/CC) 


FLOW 
PATE 
(  G  A  L  /  M  I  N  ) 


PRESSURE 
DIFFFP, 

(PSI  J 


OBSERVED 

FRICTION 
FACTOR 


COMPUTED 
FRICTION 
FACTOR 


13  7  699.     0.7761 
MAX  EXPCCTED  VALUES 

MIN  EXPECTED  VALUES 


0.7655  3,16  374.2  0.00268 
DUE  TO  UNCERTAINTY  37*. 9  0.00290 
DUE  TO  UNCERTAINTY    370.1    0.00285 


0  .0041 1 


113657.  0.7761  0.765^  2.61  223.7  0.0025-1 
MAX  EXPECTED  VALUES  DUE  TO  UNCERTAINTY  225.4  0.00255 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTAINTY    221.0    0,0025" 


0  .  00431 


86920.  0.7837  0,7655  2.03  118.7  0.00222 
MAX  EXPECTED  VAlUES  DUE  TO  UNCERTAINTY  119.8  0.00224 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTAINTY    116.0    0.00?17 


0.00461 


67440.  0.7837  0.7655  1.57  62.9  0.00195 
MAX  EXPECTED  VALUES  DUE  TO  UNCERTAINTY  63.7  0.00198 
MIN  EXPECTED  VALUES  DUE  TO  UNCERTAINTY     61.0    0.0  019" 


0.00491 


5  4  703.     0.733" 
MAX  EXPECTED  VALUES 
MIN  EXPECTED  VALUES 


0.7655  1.28  38 
HUE  TO  UNCERTAINTY  38 
DUE  TO  UNCERTAINTY     37 


0.00181 
0.00163 
0 .00175 


0.00517 
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ANALYSIS  or  flow  DATA 
FLOW  OF  50WPPM  pI-3  IN  JET  FUEL< FIRST  saSS  THRU  ORIGINAL  SYSTEM) 

in  flow  t u b e  number  i 5 . 


FLOW  TUBE  DIMENSIONS  (INCHES)   FLOW  TUBr  DIAMETER  0.0533 

TEST'SECTION  LENGTH   5.406 

REYNOLDS  COMPUTED   MEASURED    FLOW     PRESSURE  OBSERVED   COMPUTED 
NUMBER    VISCOSITY  DENSITY     RATE     DIFFPP.   FRICTION   FRICTION 
(CD)   .  (V-'/CC)        (GAL/MIM)   (PSD    FACTOR     FACTOR 


10B972.     0.736S    0.7633       1.77 

MAX  EXPECTED  VALUES  DUE  TO  UNCERTAINTY 

MlN  EXPECTEn  VALUES  DUE  TO  UNCERTAINTY 

92239.     0.7368    0.7633       1.50 

MAX  EXPECTED  VALUES  DUE  TO  UNCERTAINTY 

MIN  EXPECTED  VALUES  DUE  TO  UNCERTAINTY 

73246.     0.7368    0.7633       1  .  1 9 

MAX  EXPECTED  VALUES  DUE  TO  UNCERTAINTY 

MIN  EXPECTED  VALUES  DUE  TO  UNCERTAINTY 

56841.     0.7368    0.7633       0.Q3 

MAX  EXPECTEn  VALUES  DUE  TO  UNCERTAINTY 

MIN  EXPECTED  VALUES  DUE  TO  UNCERTAINTY 

300.     0.7366    0.76  3  3       0.00 

1362.     0.7387    0,7653       0.02 


471, 
4  7  3, 

465, 

,1 
,9 
,9 

0  .00348 
0 .00350 
0.00344 

0, 

,00435 

286, 
28^ 
283  , 

,  4 
,8 
,6 

0.00295 
0.0029Q 
0 .00292 

0  , 

,00454 

175, 

177 

172 

,  4 
.  2 
.0 

0.00287 
0 . 00290 
0.00281 

0, 

,00481 

74 

74 
72 

.  0 
.9 
.0 

0.00201 
0.00203 

0 .00196 

0, 

.00512 

o, 

.1 

0  .06749 

0, 

.05328 

0 

.3 

0 .01353 

0 

.01175 
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analysis  of  flow  data 

flow  of  jpd   with  degraded  fie  (second  pass  thru  modified  system) 

in  flow  t u b f  number  i 5 . 


FLOW  TUBE  DIMENSIONS  (INCHES)   FLOW  TUBE  DIAMETER  0.0533 

TEST  SECTION  LENGTH   5.406 


REYNOLDS  COMPUTED   MEASURED    FLOW     DRESSURE  OBSERVED   COMPUTED 
NUMBER    VISCOSITY  DENSITY     RATE     DIFFER.   FRICTION   FRICTION 
(C=»)     <  G  M  /CO)   (  G  A  L  /  M  I  N  )   (  P  S  I  )    FACTOR     FACTOR 


17181. 
14744  . 

12  0  3  0. 

8917  . 

12000. 


0.7581 
0.7581 
0 .7581 
0.7581 

0  .7583 


0.7669 
0.7669 
0  ,7669 
0 .7669 
0  .7669 


0,29  23.2  0.00655  0.00691 

0.25  18.2  0.0  0698  0.00713 

0.20  13.1  0.00751  0.00755 

0.15  7.9  0.00826  0.0  0814 

0.2  0  13.4  0.0  0  774  0.0  0  756 
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analysis  of  flow  data 

flow  of  50l-.;  ppm  p  ib  in  jet  f  u  e  l  (  f  i  r  s  t  pass  thru  modified  system) 

in  flow  tube  number  15. 


FLO  4  TUBE  DIMENSIONS  (INCHES)   FLOW  TUBE  DIAMETER  0.0533 

TEST  SECTION  LENGTH   5.406 


REYNOLDS  COMPUTED   MEASURED    FLOW     PRESSURE  OBSERVED   COMPUTED 
NUMBER    VISCOSITY  DENSITY     RATE     DIFFER,   FRICTION   FRICTION 
(CP)     (GM/CC)   (GAL/MIN)   (PSD    FACTOR     FACTOR 


18487. 

19250. 
20016. 
16160. 
13579. 


0  .  762^ 
0  .762* 
0.7624 
0.7624 

0.762- 


0  .7654 
0  .7654 
0  ,7654 
0  .7654 
0  .7654 


0,31  19.8  0.00477  0.00678 

0,32  22.0  0.00487  0.00672 

0.34  2  3.1  0.00474  0.0  0  665 

0 .27  17.2  0 .00542  0  .00702 

0.23  13.7  0.0  0  611  0.00733 
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ANALYSIS  OF  riQW  DATA 

FLOW  0 ;"  5  0  W  P  P  M  p  13  IN  JET  F  U  l:  L  (  F  1  p  S  T  PASS  THRU  MODIFIED  SYSTEM) 

IN  F'LOW  TUBE  NUMBER  15. 


FLOW  TU^E  DIMENSIONS  (INCHES)   FLOW  TUBE  DIAMETER  0.0533 

TEST  SECTION  LENGTH   5.406 


REYNOLDS  COMPUTED   MEASURED    FLOW     PRESSURE  OBSERVED   COMPUTED 
NUMBER    VISCOSITY  DENSITY     RATE     DIFFER.   FRICTION   FRICTION 
(CP)     (GM/CC)   (GAL/MIN)   (PSD    FACTOR     FACTOR 


17035. 

0  .7578 

0 

,7642 

0  .28 

17  .4 

0.00498 

0  .00602 

11477. 

0  .7578 

0 

.7642 

0  .19 

in  .i 

0  .00634 

0  .00764 

8749. 

0  .7^78 

0 

,7642 

0.15 

6.8 

0  .00736 

0.00818 

6955. 

0.7578 

0 

.  76^2 

0.12 

4  .8 

0  .00824 

0.00866 

5240. 

0  .7578 

0 

.7642 

0  .09 

2.8 

0.00851 

0.30930 

4058. 

0  .7578 

0 

.7642 

0  .07 

1.9 

0  .00943 

0.00991 

3496. 

0.7^78 

0 

.7642 

0  ,06 

1.3 

0.00900 

0.01029 

2820  . 

0  .7578 

0 

.76^2 

0  .05 

0.9 

0.00940 

0.01085 

2653. 

0.7578 

0 

.76^2 

0  .04 

C.8 

0  .00889 

0.01102 

2290. 

0  .  7^78 

0 

.7642 

0  .04 

0.5 

0.00841 

0.01143 

1401. 

0.7578 

0 

.7642 

C.02 

0  .3 

0.01231 

0.01142 

875. 

0  .  7578 

o 

.764? 

0.01 

0.2 

0  .  01954' 

0.01329 

"W 


